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ABSTRACT 


In  this  report,  which  is  the  first  part  of  a more  compre- 
hensive work,  we  pose,  and  give  a solution  algorithm  for,  a hyper- 
bolic conservation  law.  The  algorithm  is  of  an  embedding  type,  in 
that  the  solution  is  built  up  from  Green's  functions  for  simple 
processes  taking  place  without  reference  to  boundaries,  and  the 
locations  of  shocks  are  not  explicitly  followed.  There  is  a dis- 
cussion of  boundary  conditions  and  treatment  of  the  Burgers  and 
Korteweg-de  Vries  equations.  Cursory  mention  is  made  of  the 
extension  to  systems.  Appropriate  function  spaces  for  solutions 
are  Introduced.  Tlie  effects  of  perturbations  in  the  initial  con- 
ditions and  of  the  velocity  of  propagation  of  disturbances  are 
analyzed.  For  a monotonic  velocity  profile,  in  the  one-dimensional 
case  without  boundaries,  a convergence  proof  and  error  estimate 
are  given. 
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FOREWORD 


This  report  and  the  ones  to  follow  will  comprise  the  bulk 
of  a work  entitled  "Algorithms  for  Hyperbolic  and  Hydrodynamic 
Free  Boundary  Problems".  I gratefully  acknowledge  the  support 
for  this  work  provided  by  the  Office  of  Naval  Research  under  Task 
No.  NR  334-003.  This  task,  with  the  title  "Ship-Wave  Interactions", 
has  as  Its  purpose  the  computation.  In  an  efficient  and  reliable 
manner,  of  the  phenomena  attendant  to  the  Interaction  of  a rigid 
body  with  water  possessing  a free  surface.  It  Is  also  a pleasure 
for  me  to  Indicate  the  time  I spent  at  the  Applied  Mathematics 
Research  Institute,  sponsored  by  the  Office  of  Naval  Research  under 
Contract  No.  N00014-75-C-0921  with  the  Applied  Institute  of  Mathe- 
matics, Inc.,  at  Dartmouth  College  during  the  summers  of  1975  and 
1976,  when  parts  of  this  work  were  done.  It  Is  my  Intention  to 
have  the  whole  work  published  separately  under  one  cover  once 
the  various  parts  have  been  assembled. 

The  algorithms  presented  here  are  all  seml-analytlcal.  In 
that  time,  but  not  space.  Is  discretized.  Solutions  to  free  bound- 
ary problems  are  built  up  out  of  solutions  to  simple  linear  partial 
differential  equations.  The  particular  type  of  quadrature  scheme 
used  to  solve  these  simple  equations  approximately,  and  thereby 
effect  the  final  numerical  Implementation  of  the  algorithms  on  a 
computer.  Is  left  to  the  discretion  of  the  reader.  All  proofs  of 
convergence  refer  to  the  seml-analytlcal  algorithms  presented. 

In  the  work  to  follow,  there  will  be  some  unevenness  of 
rigor  and  completeness,  for  wlilch  we  offer  no  apology.  In  the 
first  part  of  the  work  we  present  some  algorithms  for  the  solution 
of  free  boundary  problems.  The  hyperbolic  problems  are  not  moti- 
vated so  much  physically  as  they  are  mathematically,  to  illustrate 
the  Ideas  that  we  bring  to  bear  on  the  hydrodynamic  problem.  There 
Is  more  physical  motivation  for  the  algorithm  for  hydrodynamic  flows, 
and  we  present  the  motivation  with  the  algorithm.  More  recondite 
features  of  the  algorithm,  however,  are  treated  more  fully  only  at 
the  end  of  this  work. 

In  analyzing  the  well-posing  of  the  hydrodynamic  Initial 
value  problem  In  the  sense  of  Its  unique  and  stable  solvability  by 
the  computational  procedure  to  be  presented,  we  will  need  to  clari- 
fy the  various  types  of  Instabilities  which  can  occur.  Among  the 
flows  which  can  evolve  In  this  way  from  the  Initial  conditions  are 
those  which  may  conventionally  be  regarded  as  "turbulent".  Ue  will 
show  the  relationship  between  stability,  energy  conservation,  and 
turbulence.  Accordingly,  a considerable  portion  of  the  work  will 
contribute  toward  the  outline  of  a rational  theory  of  Invlscld  In- 
compressible flow,  especially  the  free  boundary  problem,  and  the 
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CHAPTER  ONE 

A HYPERBOLIC  FREE  BOUNDARY  PROBLEM 


Statement  of  the  Problem  and  Presentation  of  a Solution  Algorithm 
We  start  by  considering  the  initial  value  problem 

Ut  + v(u)u^=0  , .»  < X < » , t>0,  (1.1.1a) 

U(X,0)  = Uo(x)  , -CO  < X < 00  . (1.1.1b) 


is  is  a model  hyperbolic  initial  value  problem,  and  we  will  present 
a semi-analytical  algorithm  for  its  solution.  The  proof  of  conver- 
pnce  of  the  algorithm,  for  important  classes  of  functions  v(-)  will 
be  presented  in  Chapter  Two.  ' '' 


known,  there  is  generally  no  guarantee  that  problem 
(1.1.1)  possesses  a unique  solution,  globally  in  time.  Even  if  a 
unique  solution  exists  for  an  initial  time  interval,  there  will  be  a 
tendency  for  internal  discontinuities  ("shocks")  to  develop.  After 
the  appearance  of  a "shock",  equation  (1.1.1)  can  be  solved  only  in  a 
generalized  sense,  and  must  be  supplemented  by  various  physically  moti- 
vated  8|^ck  conditions"  which  are  prescribed  at  each  such  disconti- 
nulty.  Thus,  the  problem  becomes  a free  boundary  problem,  since  the 
equation  (1.1. p will  at  best  possess  classical  solutions  only  through- 
out domains  which  do  not  contain  shocks,  but  which  have  shocks  as 
oundaries,  and  the  locations  of  the  shocks  must  be  determined  as 
part  of  the  solution  of  the  problem. 

Since  problem  (1.1.1)  is  not  in  itself  physically  motivated 
artificially  provide  an  underlying  "physical"  principle  by  re- 
casting it  in  "conservation"  form: 


‘'t  “ 0 • 

(1.1.2a) 

F(u)  - J*“  v(§)d5  . 

(1.1.2b) 

Assume  that  uodx  and  t^(x)dx 


are  defined,  and  that  u(x)  - 0 
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i; 


as  I xj  -•  ®.  We  get 


I 

I 

( 

i 


i 


Ui  =^.00  u(x,t)dx  » Uo(x)dx 


(1.1.3) 


In  developing  an  algorithm  to  solve  (1.1.1),  we  will  take  care  to 
conserve  in  (1.1.3).  Moreover,  we  shall  supplement  the  p'oblem 
by  the  condition  that 


^Ua(t) 


s 0 


(1.1. ^a) 


where 

Ug  = f (x,t)dx  . 

00 


(1.1.4b) 


(It  should  be  pointed  out  that  the  conditions  (1.1.3)  and 
(1.1.4)  serve  only  as  mathematical  motivation,  so  that  we  may  illus- 
trate, for  the  simple  problem  (1.1.1),  ideas  which  will  be  developed 
more  fully  later.  In  fact,  we  will  see  in  a later  report  that  an 
equation  like  (1.1.1a),  with  v(u)  - u,  arises  in  a very  simple  flow 
problem,  but  that  the  equation  by  itself  is  in  the  most  general  case 
devoid  of  meaning  unless  it  is  coupled  to  an  equation  + (pu)  <=  0. 
In  this  case,  the  underlying  physics  of  the  problem  shows  that  the 
correct  generalized  solution  of  (1.1.1a)  depends  critically  on  the 
function  p(x,0)  = po(x).) 

The  algorithm  we  use  to  solve  (1.1.1),  supplemented  by  (1.1.3) 
and  (1.1.4),  is  of  an  "embedding"  type,  in  which  the  shocks  are  not 
followed  explicitly.  Such  a method  is  very  efficient  to  treat  cases 
where  there  are  many  shocks,  or  when  there  are  several  independent 
"space"  variables  x,  as  we  shall  see  in  the  following.  The  technique 
has  something  in  comnon  with  other  embedding  techniques  like  the 
method  of  artificial  viscosity  (Reference  11),  except  that  here  no 
extraneous  parameter  is  introduced,  beyond  the  time  step  t in  terms 
of  which  the  evolutionary  problem  is  solved.  The  method  we  use  is 
built  on  the  basic  conservation  "law"  (1.1.3),  which  we  regard  as 
more  fundamental  than  the  differential  equation  (1.1.1a).  Similar 
embedding  methods  based  on  conservation  laws  have  been  successfully 
used  for  other  free  boundary  problems  (References  2 and  3). 

Before  proceeding  with  an  algorithm  to  solve  (1.1.1),  let 
us  consider  a slightly  generalized  version  of  the  problem,  where  the 
solution  of 


*1 

f 


“..■••(FCu))  "0,  Xo<x<Xi,  t>0,  (1.1.5a) 

c X 
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u(x,0)»  Uo(x),  Xo  S X s Xi  , (1.1.5b) 

supplemented  by  boundary  conditions  at  Xq  and  x^  , is  sought.  (It 
turns  out  that  some  natural  boundary  conditions  for  the  problem  are 
expressed  simply  in  terms  of  a distribution  function  which  appears 
in  the  execution  of  the  solution  algorithm.  Therefore,  an  explicit 
statement  of  the  boundary  conditions  is  deferred  until  we  describe 
the  solution  algorithm  in  detail,  below.)  We  thus  have  an  initial- 
boundary value  problem. 

Generalizations  of  equations  (1.1.3)  and  (1.1.4)  are 


dt 


u(x,t)dx  = -F(u(b,t))  + F(u(a,t)), 


Xo  ^ a < b ^ Xj^  , 


and 


(1.1.6) 


u2(x,t)dx  £ - I^(u(b,t))  + Fs  (u(a,t)>. 


XqS  a<  b < Xi  , 


where 


(1.1.7a) 


Fs(u)  = 2 5 v(?)d? 


(1.1.7b) 


Before  we  proceed  with  a precise  description  of  our  solution 
algorithm  for  the  Initial  and  boundary  value  problem  associated  with 
(1.1.5),  let  us  give  a schematic  description.  To  facilitate  this 
description,  we  shall  assume  that  u(x,t)  a 0 for  all  x € [xo,x^].  We 
may  assume  we  have  a graph  of  the  function  u(-,t),  for  t fixed,  as 
shown  in  Figure  1. 


Figure  1.  Schematic  Solution  u (•.  t),  for  t Fixed. 
-11- 
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Since  in  our  conservation  law  (1.1.6)  changes  in  the  quan- 

-b 

tity  u(x,t)dx  are  brought  about  through  terms  which  are  functions 

only  of  u at  the  end  points  a and  b,  it  is  suggestive  for  us  to  think 
pb 

of  J u dx  as  the  total  amount  of  some  "fluid"  which  is  conserved, 
a 

except  for  the  amount  "flowing"  across  the  vertical  lines  x = a and 
X = b.  In  order  to  avoid  confusion  of  notation,  we  think  of  the 
plane  in  Figure  1 as  the  xo-,  instead  of  xu-,  plane.  The  quantity 


p b 

J u(x,t)dx  is  just  the  area  in  the  xo-plane  contained  between  the 

curves  o'  = 0 and  O'  = u(x,t)  and  the  vertical  lines  < = a and  x = b. 
Thus,  we  may  think  of  the  flow  as  area-preserving,  except  for  amounts 
crossing  x = a and  x = b.  One  of  the  principal  fe«>tures  of  our  al- 
gorithm is  that  this  area  is  conserved,  except  for  such  flow  across 
the  boundary. 


In  addition,  we  may  think  of  the  "fluid"  as  comprised  of 
"streams",  each  moving  along  at  constant  a.  The  portion  of  fluid 
contained  between  a and  a + da  and  the  vertical  lines  x = a and 
X = b has  the  "area" 


where 


pb 

da  ^ f (x, t,a)dx 
a 


f(x,t,a) 


1 OSaiu(x,t) 
0 a > u(x, t) 


(1.1.8) 


The  whole  fluid  area  between  x = a and  x ■>  b is 


“O  •'a  • 

At  each  point 

u(x,t)  - j f(x,t,a)da  . (1.1.9) 

0 

Hence  the  fluid  flow  may  be  envisaged  as  the  flow  of  an  In- 
finite number  of  "streams",  labeled  by  the  parameter  a,  each  with 
its  own  area.  The  area  of  each  stream  will  change  due  to  the  flow 
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of  fluid  in  Che  stream  across  Che  spatial  boundaries,  and  also 
due  to  transfer  of  fluid  from  one  stream  to  another.  The  flow 
of  fluid  in  a stream  across  the  spatial  boundaries  will  be 
governed  by  the  boundary  conditions.  The  transfer  of  fluid  from 
one  stream  to  another  will  take  place  in  such  a way  that  the  total 
area  in  the  streams  is  conserved.  The  rules  governing  the  manner 
in  which  fluid  is  transferred  from  one  stream  to  another  inter- 
nally will  constitute  an  expression  of  the  shock  conditions  which 
are  imposed.  The  various  streams  are  indicated  in  Figure  1 by 
Che  horizontal  lines. 

To  implement  our  solution  algorithm  we  first  convert  each 
stream  independently  according  to  the  linear  equation 


fj.(x,t,Qf)  + vCo-)  f^(x,t,a)  = 0 


(1.1.10) 


for  a small  time  interval,  called  a time  step;  Equation  (1.1.10) 
can  be  solved  exactly.  We  sweep  through  the  streams  in  the  order 
of  increasing  or. 

When  v(o)  > 0 (v(a)  < 0)  the  stream  between  a and  a + do 
enters  (exits)  the  flow  domain  at  xq  and  exits  (enters)  at  x^ . As 
is  well  known,  for  equations  like  (1.1.10)  we  should  thus  specify 
f(x,t,o)  at  Xo  when  v(a)  > 0 and  at  x^  when  v(of)  < 0.  This  informa- 
tion will  be  provided  by  the  boundary  conditions.  When  the  range 
of  values  of  u in  the  problem  is  such  that  we  will  have  v(u)  > 0 
in  one  region  and  v(u)  < 0 in  another,  we  may  generally  expect  that 
it  will  be  necessary  to  prescribe  some  sort  of  boundary  data  at 
both  xq  and  Xj^  , but  in  such  a way  as  not  to  overdetermine  the  prob- 
lem. 


Figure  2 shows  schematically  Che  domain  in  Che  xo-plane 
occupied  by  the  fluid  in  Figure  1 after  convection.  The  arrows 
point  to  the  right  for  v(cr)  > 0,  to  the  left  for  v(or)  < 0. 

If,  after  convection,  Che  streeun  corresponding  to  o » cyo 
has  flowed  into  a spatial  dcmaln  which  has  not  been  filled  with 
all  the  streams  for  a < cxq,  we  let  the  "fluid"  In  the  stream  a » 
Qfo  "fall"  or  "cascade"  down  to  occupy  the  smallest  unfilled  values 
of  Of  < «o  each  point  in  the  domain,  in  such  a way  that  the 
total  area  of  the  streams  remains  unchanged.  The  result  is  shown 
in  Figure  3. 

A more  mathematical  description  of  our  method  for  obtain- 
ing an  approximate  solution  to  the  problem  (1.1,5),  subject  to 
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boundary  conditions  to  be  specified,  follows.  We  shall  eliminate 
the  restriction  u(x,t)  ^ 0 made  In  our  heuristic  description  above. 
In  addition  to  the  Initial  function  Uq(x),  we  require  functions 
u(xo>t)  and  u(Xi,t)  to  be  prescribed.  These  functions  will  be 
used  to  supply  the  boundary  data  required  to  solve  Equation  (1.1.10) 
uniquely,  as  described  above.  We  assume  that  Uq(x),  u(xo,t),  and 
u(xi  ,t)  are  bounded  from  below: 


1 


u(x^,t)  i u<,  1=0,  1,  t € (0,Tl  , 

Uo(x)  i u^,  X € [xq  ,xi  ] , 


(1.1.11) 


for  some  real  number  u^.  Let  t > 0 be  given.  We  start  with 

" Uo(x),  x6  [xo,Xil  , (1.1.12) 


and  the  algorithm  will  be  complete  when  we  specify  how  to  get  from 
u*'(x)  to  u^"^i(x).  Consider  the  solution  of  the  linear  problem 


(x.t.Qf)  + v(Qf)  f"  (x,t,a)  = 0 (1.1.13a) 

subject  to  the  Initial  condition 


f"(X,0,CK) 


0 or  < u< 

1 u<  sq;s:u"(x),  x€[xo.Xi1  (1.1.13b) 

0 a > u"(x) 


and  the  boundary  conditions 


f"(Xo,t,o)  - 


0 

1 

0 


a < u^ 

s asu(j9,,t+ni)  , v(a)  >0,  0 < t , 

a > u(xo,t+nT),  v(a)  > 0 


(1.1.13c) 
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0 a < 

1 ^ a £ u(xi  , t+tiT)  , v(a)  <0,  0<t.  (1.1. 13d) 

0 a > u(xj  jt+nr),  v(o')  < 0 

Finally,  let 


f”(Xi  ,t,Qf)  = 


u"^^(x)  = U + J f"(X,T,Q')dO'  , x€  [Xo.Xi]  . (1.1.14) 

^ - 00 


It  is  our  hope  that  in  some  sense  u*'(x)  will  approximate  the  exact 
solution  u(x,nT)  of  problem  (1.1.5),  supplemented  by  suitable  bound- 
ary conditions.  We  call  the  parameter  t the  time  step.  In  the  next 
chapter  we  will  prove  the  convergence  of  u*'(x)  to  u(x,t),  as  t =• 
t/n  i 0,  or  n t “,  for  Important  classes  of  problems  of  the  sort 
(1.1.5),  and  we  will  obtain  a bound  on  the  error  involved  in  approxi- 
mating u(x,t)  by  u"(x). 

Note  that  the  equation  (1.1.13)  for  f”  may  be  solved  explicitly. 
The  solution  is,  for  0 £ t. 


f"(x,t,o)  . 
where 


f‘'(x,0,cr),  x6  Cxo,Xj],  or  ^ Cq  , (1.1.15) 

f (max(xo,x-v(a)t),  max(0,  t-  , x€[xo,x^J,  a € , 

f"(min(xi  ,x-v(o')t  ),  «nax(0,  t-  [xq  ,Xi  ],  or  € C_ 


Co-  (or|v(a)  -0},  {olv(a)>0},  C_-  {«!  v(o)  < 0}  , (1.1.16) 


Thus,  at  no  step  in  the  implementation  of  the  algorithm  need  we  take 
finite  differences.  The  expression  (1.1.15)  may  be  thought  of  as 
giving  the  solution  in  terms  of  the  initial  and  boundary  data  by 
means  of  Green's  function  for  equation  (1.1.13a),  and  Green's  func- 
tion fox  this  problem  is  everywhere  non-negative.  Ihe  entire  algo- 
rithm is  executed  as  a series  of  integrations  of  non-negative  quan- 
tities. 

At  this  point  let  us  state  the  boundary  conditions  we  are 
imposing  on  the  problem  (1.1.5).  We  let  the  "density"  u(x,t)  be 
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related  to  a ''distribution  function"  f(x,t,a)  by 


u(x,t)  = If(x,t,‘)  = u<  + f(x,t,0')dQ',  x€[xb,Xjl,  t e [0,t]  , 

(1.1.17a) 

iO  a < u< 

1 i O'  s U(x,t),  X e [xo.x^],  t e [0,Tl 

0 a > u(x,t)  (1.1.17b) 


For  any  (x,t)  € [xo.t^]  X [0,T]  we  can  decompose  u(x,t)  as  follows: 


u(x,t)  = u(x,t)  + u (x,t)  + u (x,t) 


(1.1.18) 


where 


u(x,t)  - + j^Xo(Qf)  f(x,t,a')dor 

u (x.t)  “ Xoo  f(x,t,a)daf  , 


u (x.t)  f(x,t,or)do'  , 


(1.1.19) 


Xo(«) 


1 v(cr)  = 0 

0 v(o()  ^ 0 


X^(a) 


1 v(a)  > 0 

0 v(of)  s 0 


(1.1.20) 


X_(Qf) 


1 v(a)  < 0 

0 v(a)  i 0 


Then  the  boundary  conditions,  which  appear  in  Che  algorithm  as 
(1.1.13c,d),  in  terms  of  the  specified  functions  u(xo.t),  u(xj,t), 
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when  taken  in  conjunction  with  the  Initial  conditions  and  the  equa- 
tion of  evolution,  amount  to  specifying  lira  u'*’(x,t)  and  lira  u (x,t), 

X I Xq  X tXj 

Only  in  special  cases,  such  as  when  v(or)  > 0 for  all  oi,  so  that 
X^(q()  = 1,  do  we  find  that  this  corresponds  to  prescribing  u at  one 
or  the  other  of  the  end  points. 

A relation  like  (1.1.17a)  may  be  used  in  conjunction  with 
the  expression  (1.1.15)  for  f”(x,t,cr)  to  define  a function  u”(x,t): 


u"(x,t)  - I (f"(x,t-nT, •)) 


(1.1.21) 


Note  that,  according  to  (1.1. 15) , (1. 1. 13) , and  (1.1.14), 


u"(x,nT)  » u"(x)  , 


X 6 [ Xo  , Xi  1 , 


(1.1.22) 


u"(x,(n+l)T)  - u"'*’^(x),  x€  [xo,Xj]  . 


With  u'^(x,t)  defined  in  this  manner,  we  can  show  that 
results  similar  to  (1.1.6)  and  (1.1.7)  are  obtained,  for  t-nr 
sufficiently  small  and  suitable  requirements  of  regularity  on  the 
initial  function  Uq(x)  and  the  boundary  functions  u(x^,t).  To 
show  this,  we  first  note,  from  equation  (1. 1. 15),  that*^ 


f”(x,t,Qf)dx  f"(x,0,Of)dx  + v(of)  f"(a,t' ,Q;)dt' 


rb  ,n 


(1.1.23) 


v(of)  f"(b,t' ,Of)dt' 


Integrate  this  over  or  and  use  (1.1.21),  (1.1.13b),  and  (1.1.2b): 


u"(x,t+nT)dx-^**  u”(x)dx  + t F(u"(a))-  t F(u"(b)) 


(1.1.24) 


■‘•4^Hv(or)  [f"(a,t',o)-f"(a.O,of)-f"(b,t',o)+f"(b.O,o)]<to  dt' 
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Differentiating  (1.1.24)  with  respect  to  t,  we  get 


^ j’*’  u"(x,t+nT)dx  ■ F(u“(a))  - F(u"(b)  ) 


+ j*“  v(or)  [f"(a,t,a)-f"(a,0,a)-f"(b,t,cr)+f"(b.0.o)lda 


In  order  to  bound  the  last  expression  on  the  right  of 
(1.1.25),  we  will  find  it  convenient  to  introduce  some  more  nota- 
tion. Let 


U = max 


/ sup  Uo(x),  sup  u(xo,t),  sup  u(xi  ,t)\  (1.1.26) 

\x€[xo,Xil  t€(0,T]  t€(0,T]  / 


From  equations  (1.1.12)  - (1.1.14)  and  (1.1.21),  we  see  that 


u"(x,t)  s U and  f"(x,t,o()  - 0 for  a > U . (1.1.27) 


Furthei;  define 


V = sup  lv(a)l 
a€[u^,U] 


(1.1.28) 


Consider  a representative  contribution  to  the  right  hand  side  of 
(1.1.25): 


X*  v(o)  [f"(a,t,Of)  - f”(a,0,o)]dof  . 


If  xo  & a - v(of)t  s Xj  , write 


•f"(A|t.»)  - f"(a,0,of)  - f“(a-v(o)t,0,a)  - f"(a,0,o)  (1.1.29) 


YHE  X)HNS  »0PKINS  university 

applied  physics  laboratory 

LAuRfl 


by  (1.1.15),  and  note,  from  (1.1.13b)  that  the  right  hand  side 
of  (1.1.29)  vanishes  unless  a lies  between  u*^(a)  and  u”(a-v(or)t) . 
Thus,  suppose  Xq  + Vt  s a s - Vt.  It  follows  from  (1.1.27) 
and  (1.1.28)  that  Xq  s a - v(o)t  ^ x^  for  all  a such  that  the 
difference  in  (1.1.29)  is  non-zero.  We  find 


P"  U" 


lf"(a,t,or)-f"(a,0,a)|do  s sup  (u"(§i  )-u"(52) ) 

?i  ,§a€[a-Vt,a-+Vt] 


in  this  case.  If  a < Xq  + v(a)t,  write 


(1.1.30) 


f"(a,t,a)  - f"(a,0,o)  - f"(xo,0,a)  - f"(a,0,a) 


^ (xo>t“  v(cr)  ^ (^>0,Qf)  . 


(1.1.31) 


By  the  same  reasoning  that  led  to  (1.1.30),  we  get,  for  a < Xq + Vt 
and  a s x^  “ Vt, 


j |f”(a,t,Qr)-f"(a,0,or|da  s sup  (u"(5i  )-u"(52) ) 

.5a€[xo,a+Vt] 


(1.1.32) 


(’^)"U^(*b.Tl+nT))  , 

where  (u  ) and  u are  obtained  in  terms  of  u*'  and  u through  rela- 
tions like  (1. 1. 17)-(1. 1.20) . The  remaining  contributions  to  the 
right  hand  side  of  (1.1.25)  are  bounded  similarly,  and  we  get 


jJ*  u(x,t+nT)dx  - F(u"(a))  + F(u”(b))  | s:  vo  (t) 


(1.1.33a) 
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where 

On(t) 


+ 


+ 


+ 


5i  •?a€[max(xo  ,a-vE)  ,min(xi  ,a+Vt)]  (5a)) 

5i  ,5a6[max(xo  ,b-V?)  ,min(Xj  ,b-Wt)]  ^ “ (5a)) 


(1.1.33b) 

TiI[8  tl  + ^8jP^^((«V(x<i>-5'*(xb.T)fnT)  ) 

+ ^g®5Pj.^((“")’(’‘i)-2‘(xi  ,11+nT))  . 


An  Integrated  version  of  this  equation  becomes,  on  use  of 

(1.1.22), 

1 u (x)dx  - Uo(x)dx  - T 2^  [F(u^(a))  - F(u^(b))]l  5 Vt ^2^  0^(7)  . 

(1.1.34) 


The  right  hand  side  of  (1.1.33)  will  go  to  zero  as  x 1 0 
for  0 < t s T , if  we  have 

(i)  the  u are  continuous  in  some  neighborhoods  of  a and  b; 

(ii)  the  boundary  functions  u'^(xo,‘)  and  u"(xi,‘)  are  contin- 
uous in  a neighborhood  of  nx  + t; 

(ill)  (u  ) (xo)  ■ lira  u^(xo,t)  and  (u*')  (x.)  - 11m  u (Xj  ,t)  . 

tlnx  tinx 

Those  regularity  properties  which  relate  to  the  calculated  functions 
u will  often  follow  from  an  a^  priori  or  a posteriori  analysis  of 
the  problem.  For  example,  in  the  important  special  case  v(o)  > 0, 
o € Cu^,U],  we  have  (ill)  above  for  n > 0.  In  the  case  a • xq  "* 

-®>  and  b ■ x^  ->+<*>,  for  most  types  of  asymptotic  conditions  of 
interest,  (i)  holds.  There  will  be  mention  of  generalized  ver- 
sions of  (1.1.33),  for  different  classes  of  initial  and  boundary 
data,  in  Chapter  Two. 
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A result  similar  to  (1.1.7)  follows  if  we  multiply  (1.1.23) 
by  a and  integrate  over  a.  Note  that 


ia.  “ f"(x,0,Qf)da  [(u”(x))^-u^  j (1.1.35a) 

and 

la.  “ ^ 7 [("'*(*»  • (1.1.35b) 

Thus 

jj’  (u"(x,t+nT)j2dx  - (u"(x))®dx  - tF2^u"(a))  + tFs  (u"(b))  S G 

(1.1.36a) 

where 


G = 2 J_^a  v(0f)  [f”(a,t',Of)  - f"(a,0,a) 


(1.1.36b) 


- f"(b,t',a)  + f"(b,0,cr)]dadt'  . 


Proceeding  as  above,  we  get 


1g1  s 2t(|u|  + lu^l)  Vo^(t) 


(1.1.37) 


The  analogue  to  (1.1.34)  is  obtained  by  using  (1.1.22)  and 
(1.1.37): 

(u"(x))®dx  - (uo(x))=dx  - T |Fa  (u^(a))  - Fa  (u^(b))j  s H 

(1.1.38a) 

where 

. n-1 

|h|  s 2t(1u1  + lu^l)  V ^ o^(T)  . (1.1.38b) 
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The  right  hand  side  of  (1.1.37)  will  be  o(t)  as  t i 0, 
under  conditions  (i)  - (iii)  given  above. 

In  view  of  the  relations  (1.1.33),  (1.1.34),  (1.1.36), 
(1.1.37),  and  (1.1.38),  we  may  think  of  the  algorithm  (1.1.12)- 
(1.1.14)  as  "conserving"  the  quantity  given  in  (1.1.3),  and 
"diminishing"  Ug  given  in  (l.I.4b),  when  those  quantities  are 
defined.  It  is  natural,  then,  since  and  are  given  as 
integrals,  to  Introduce  the  following  set  functions  for  Lebesgue- 
measurable  sets  (1  in  xor  space: 


Ui(f;n)  = f(x,o)  da  dx  , (1.1.39a) 

Us(f;n)  s 2 |^af(x,a)da  dx  , (1.1.39b) 


(Note  that  these  definitions  depend  on  the  value  of  the  number  u 
given  in  (1.1.11).) 

We  should  point  out  that,  if 


(i) 

(ii) 

(iii) 

(iv) 

(V) 


uq(x)  is  continuous  for  x € [xo,x^]  , 

u(xo,t)  is  continuous  for  t £ (0,5ol 
Uo(Xo), 

u(x^,t)  is  continuous  for  t £ (O,5o3  and 
Uo(Xl)» 


lim  u(xo,t)  = 
tiO 


lim  u(x^ , t)  = 
tiO 


X + v(uq(x))§  is  nondecreasing  in  ' 

X € [xo.Xjl  for  Os?  s 5o  . 

max[v(u(xo , t))  (?-t),0]  is  nonincreasing  in  t € [0,lol 
for  t s ? s ?o, 


(vi)  min[v(u(xi , t)) (?-t) ,0]  is  nondecreasing  in  t € [0,5q] 

for  t s 5 s ?o. 


then  for  0 s nr  s 5o»  the  quantity  u (x)  computed  by  (1.1.12)- 
(1.1.14)  will  be  the  exact  solution  at  time  t ■ nr  of  (1.1.5)  with 
the  boundary  conditions  u+(xoft)  ■ u+(xo,t)  and  u"(x^,t)  - u”(xi,t), 
The  quantity  u"(x,t)  given  in  (1.1.21)  will  be  the  exact  solution 
at  time  t if  t s ?g. 
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To  satisfy  ourselves  as  to  the  validity  of  this  assertion, 
we  shall  outline  a proof  for  the  case  when  t * t s ^ ““> 
Xi  — +“>  , so  that  we  only  need  conditions  (i)  and  (iv)  above. 

The  exact  solution  of  (1.1.1)  at  time  t = t is  given  by 


u(x,t)  = Qf(x)  = suptala-uo  (x-v(a)T)  s O}.  (1.1.41) 


To  show  that  u^(x)  given  by  ( 1. 1. 12) - (1 . 1. 14)  is  Qf(x) , we  will  make 
use  of  the  following  lemma. 

Lemma  1.1.1:  If  a ^ u^,  f°(yj^,0,a)  = 0,  and  f°(y2,0,a)  = 1, 

then  there  exists  y*  ^ [«nin(yi  ,ya)  ,max(yj  ,yg)]  such  that  a = u©  (y*) . 
Also,  y*  ^ y^  . 


Proof:  The  result  follows  from  (1.1.12),  (1.1.13b),  and  the 

continuity  of  u© , since  a > Uo(yi)  and  a ^ Uo(y8). 

To  show  that  u^(x)  is  Qf(x) , we  show  first  that  f°  (x-v(of)T  ,0,a) 
• 1 for  u^  s o < a(x).  We  do  this  by  contradiction.  Suppose  there 
is  OTo.  ^ Qfo  < Q'(x),  with  f°  (x-v(o(o)t  , 0,010)  = 0.  Then  Oo> 

Uq  (x-v(ofo)T) . By  the  definition  of  or(x)  , equation  (1.1.41),  there  is 
> Oq  ^ Qf(x) , with  - uq  (x-v(ofi  )t)  s 0.  We  have  the  follow- 

ing relations  imned lately: 

f®  (x-v(Ofo)T,0,Q'o)  = f°  (x-v(q'o)t,0,q'i  ) = 0 , 

(1.1.42) 

(x-v(0fi  )T,0,Qfj  ) « f®  (x-v(Oi  )T,0,Qfo)  “ 1 • 


We  cannot  have  v(orc)  » v(ofi ) without  violating  the  continuity  of  uo . 
Consider  the  case  v(oro)  > v(Qfj^).  Then  by  the  lemma  and  (1.1.42)3  y* 
such  that 


X - v(Oo)t  < y * s X - v(q'i)t  (1.1.43) 

and 

Uo(y*)  ■ Qfi  A (1.1.44) 
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Since  Uo(y*)  = Ox,  f°(y*,0,ai)  = 1 and 


f°(y*.0,ao)  = 1 


By  the  lemma,  (1.1.42),  and  (1.1.45),  3 y©  with 


X - v(ao)T  < y*  ^ y* 


Uo(yo)  = Ofo  . 


By  (1.1.43),  (1.1.44),  (1.1.46),  and  (1.1.47), 


y*  + v(uo(y*))T  > X s y?"  + v(uo(yi))T  , 


which,  with  y*  s y*,  violates  condition  (iv)  of  (1.1.40).  The 
case  v(q(o)  < v(Oj^ ) is  treated  similarly.  Thus,  f°  (x-v(O()t,0,Q')  = 1 
for  ^ a < ct(x).  Finally,  suppose  f°  (x-v(a)T,0,O')  = 1 for  some 
Of  > cr(x).  Then  u^  or  s Uo(x-v(a)T)  and  by  the  definition  of 
a(x),  or  s a(x),  givlne  a contradiction.  Putting  these  results  to- 
gether, we  see  that  u^(x)  » a(x)  follows  from  (1.1.14). 

One  may  think  of  (1.1.13a)  as  a "Boltzmann  equation", 
because  of  the  analogy  to  the  familiar  transport  equation  of 
kinetic  theory.  The  equation  we  consider  is  especially  simple, 
and  does  not  contain  the  nonlinear  "collision  terms"  which  gener- 
ally appear  in  the  Boltzmann  equation.  To  say  that  our  algorithm 
thus  represents  the  motion  of  a "collisionless"  fluid,  however, 
would  be  inaccurate.  Indeed,  the  transformation  which  takes  place 
at  each  time  step,  when  the  "distribution  function"  f”(x,T,or)  is 
converted  into  f"'^^(x,0,Qf)  by  the  processes  of  equations  (1.1.14) 
and  (1.1.13b),  represents  collisions  of  a rather  extreme  form. 

If  we  are  to  regard  the  transport  equation  (1.1.13a)  as  express- 
ing the  flow  of  independent  "streams"  at  different  "density  levels" 
a,  we  see  that  at  each  time  step,  by  some  undefined  internal 
"interaction"  mechanism,  the  various  streams  interact  in  such  a 
way  that  they  are  forced  into  the  lowest  "unfilled"  density  levels 
(above  u^).  Accordingly,  the  algorithm  (1. 1. 12)-( 1. 1. 14)  gives 
the  evolution  of  the  system  in  terms  of  a sequence  of  propagations 
of  non-interacting  streams,  and  intense  interactions  of  these 
s treams . 
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The  fact  that  equation  (1.1.5)  for  the  density  u may  be 
related  to  an  underlying  transport  equation  for  the  distribution 
function  f is  not  surprising,  if  one  recalls  that  the  equations 
of  fluid  mechanics  may  in  some  sense  be  derived  from  the  Boltzmann 
equation  (Reference  5). 


2. 


Higher-Dimensional  Cases 

A higher-dimensional  analogue 

Uj.  + V • (F(u))  = 0, 

u(x,0)  = Uo(x) , 


of  (1.1.5)  is 
X € 3),  t > 0, 


X € ® , 


(1.2.1a) 

(1.2.1b) 


with  boundary  conditions  to  be  specified.  By  ® we  mean  the  closure 
of  S).  Here  F(u)  is  a vector.  Its  derivative  with  respect  to  u is 
another  vector: 


v(o)  - ^ F(u)  . 


(1.2.2) 


As  to  boundary  conditions,  let  functions  u(x,t)  be  given 
for  x E S®  and  t € (0,T],  with 


u(x,t)  i:  u^,  X € B®,  t € (0,Tl  , 

Uq(x)  2 u^,  X € ® . 


(1.2.3) 


If  we  write 


OD 

u(x,t)  H If(x,t,*)  s u^  + J|^f(x,t,Qf)da, 


X e ®,  t€  [0,T],  (1.2.4) 


f(x,t,a)sD(u(x,t) ;a)  h 


0 or  < u^ 

1 u^  < or  s u(x,t),  x€3),  t€[0,T], 


(1.2.5) 


1 0 a > u(x,  t) 
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and  define 


C(x,t)  = {o'lx-v(a)5  € ® for  all  5 € x € S , 

C'(x,t)  H {alo  k C(x,t)],  X e S)  , 


we  can  decompose  u(x,t)  as 


u(x,t)  “ u(x,t,5)  + u'(x,t,5) 

where 

00 

u(x,t,5)  = U +J  x(x,|,0f)  f(x,t,Qf)da  , 

^ as  CO 


u'(x,t,§)  =»  x'(x,§.a)  f(x,t,Qf)dof  , 

and 

t 1 a € C(x,|) 

X(x,§,Qf)  = < 

( 0 a i C(x,5) 

x'(x,5,of)  » 1 - X(x,|,a)  . 

For  t G (0,t1,  X € 9®,  and  Of  € C^(x,|)  for  all  5 > 0,  we  set 


f(x,t,or)  - DCuCx.t);^)  . ( 


The  boundary  conditions  (1.2.10),  when  combined  with  (1.2.1), 
to  specifying 


inf  lim  u'(y,t,5),  x€  9®,  t€  (0,Tl  . ( 

§>0  y€®-x 


To  solve  (1.2.1)  subject  to  the  boundary  conditions 

we  set 


(1.2.6) 

(1.2.7) 

(1.2.8) 

(1.2.9) 

1.2.10) 

amount 

1.2.11) 

1.2.11), 
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u°(x)  = Uo(x),  X € X , (1.2.12) 

and  show  how  to  get  from  u*'(x)  to  u”"*''(x).  Given  u'^(x),  construct 

f"(x,0,Of)  = D(u"(x);q'),  X € X (1.2.13a) 

and  solve 

f"  + v(q')  • 7 f"(x,t,a)  = 0 (1.2.13b) 

subject  to  the  initial  condition  (1.2.13a)  and  the  boundary  condi- 
tion 

f"(x,t,Q')  = f(x,t+nT,Qf),  X € 3X,  t > 0,  a € c'(x,|)  for  all  ? > 0, 

(1.2.13c) 

where  f (x,  t-HiT ,of)  is  given  by  (1.2.10).  Finally,  set 

u”"'’^(x)  » If"(x,T,-),  x€  X . (1.2.14) 

The  problem  (1.2.13)  for  f may  be  solved  explicitly; 

f"(x,t,of)  - f"(x-v(Qr)[t-max(0,inf[§lci'€  C(x,t-5)})], 

, , (1.2.15) 

max(0,inf{§|Qfe  C(x,t-5)}),a)  , 

X € X,  t 2 0.  As  in  the  one -dimensional  case,  the  exact  solution 
(1.2.15)  may  be  approximated  nianerically  by  means  of  quadratures 
involving  only  non-negative  quantities. 

3.  Burgers  and  Korteweg-de  Vries  Equations 

The  method  described  may  be  used  to  solve  more  general 
equations,  and  we  indicate  how  one  would  proceed  with  some  repre- 
sentative examples.  The  extension  to  still  other  types  of  equa- 
tions will  then  be  obvious. 
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Let  us  consider  "Burgers'  equation" 


(Reference  4) 


with 


and 


+ 7 • (F(u))  =v4u,  x€®,  t>0 

u(x,0)  « uo(x) , X € Si  , 
u(x,t)  * u(x,t),  X € a®,  t € (O.Tl  . 


(1.3.1a) 

(1.3.1b) 

(1.3.1c) 


(A  more  physically  motivated  generalization  of  Burgers'  equation 
to  a system  of  equations  for  a hlgher-dlmenslonul  flow,  with  the 
flow  driven  by  the  requirements  of  conservation  of  mass  and  momen- 
tum, will  be  treated  In  a later  report  on  this  work.) 

As  above,  let  u^  satisfy 


s uo  (x)  , X € S , 


(1.3.2) 

« u(x,t),  X € a®,  t e (o,t]  . ’ 

iP  (x)  Is  given  by  (1.2.12).  Let  f''(x,t,cy)  satisfy  the  equation  (v 
given  by  (1.2.2)) 


f"  + v(0f)  • 7 f"  - V A f",  X € ®,  t > 0 (1.3.3a) 

with  f (x,0,a)  given  by  (1.2.13a)  and  the  boundary  condition 

f"(x,t,o)  - D(u(x,nT+t);o),  x € a®,  t > 0.  (1.3.3b) 

D(u;a)  Is  given  by  (1.2.5).  u"^^(x)  Is  given  In  terms  of  f*'(x,T,of) 

by  (1.2.14). 

Note  that  we  may  combine  (1.2.14)  with  the  solution  of 
equation  (1.3.3a)  subject  to  the  Initial  condition  (1.2.13a)  and 
boundary  condition  (1.3.3b)  approximately  as  follows: 
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h(x,T)  (1.3.4) 

where 

*vAh,  x€D,  t>0, 

h(x,0)  - Ig(x,T,*).  X € i)  , (1.3.5) 

h(x,t)  - u(x,nT+t),  X € 9S),  t > 0 , 

and 


g(.  + v(cr)  • V g = 0,  X € D,  t > 0, 

8(x,0,a)  - f‘'(x,0,a)  , (1.3.6) 

?<x,t,a)  « D(u(x,nT+t);or),  x€  9D,  a€c'(x,|)  for  all  | > 0,  t > 0 


The  device  of  approximately  solving  the  initial-boundary  value 
problem  for  f"  and  Integrating  to  get  by  solving  (1.3.6)  and 

(1.3.5)  is  analogous  to  the  method  of  "fractional  time  steps".  If 
2D  is  all  space,  this  method  becomes  exact,  because  of  the  commuta- 
tivity of  the  operators  vA  and  v(a)  • V.  (Solution  of  either 
initial  value  problem  (1.3.6)  or  (1.3.5)  is  then  equivalent  to 
multiplying  the  spatial  Fourier  transform  of  the  dependent  func- 
tion by  an  algebraic  factor.)  Thus,  if  we  replace,  in  the  algo- 
rithm (1.2.12),  (1.3.3a),  (1.2.13a),  (1.3.3b),  and  (1.2.14),  equa- 
tions (1.3.3a),  (1.3.3b),  and  (1.2.14)  by  (1.3.6),  (1.3.5),  and 
(1.3.4),  so  that  the  modified  algorithm  is  described  by  (1.2.12), 
(1.2.13a),  (1.3.6),  (1.3.5),  and  (1.3.4),  the  solution  of  (1.3.1) 
is  seen  to  be  effected  through  a sequence  of  convections  of  non- 
interacting  streams.  Interactions  of  these  streams,  and  diffusions. 
This  is  similar  to  the  description  given  at  the  end  of  section  1.1 
of  the  algorithm  ( 1. 1. l2)-( 1. 1. 14) , except  that  now  a diffusive 
process  is  interspersed  with  the  others. 

Consider  the  "Korteweg-de  Vries  equation"  (Reference  1) 


+ V • (F(u))  + PAu^-0,  P>0,  x€J),  t>0,  (1.3.7a) 

with 

u(x,0)  - Uo(x),  x€»,  (1.3.7b) 
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and 


where 


- S(x.t). 

-,<x.t)  . s 


X € 

X € as)"^, 

^ • V > 


‘ ^ (o.r]  , 

^ ^ (0,T1  , 

0} 


and 


n(§)  la  the  unit 


(1.3.7c) 
(1.3. 7d) 

(1.3.8) 


"«">«  to  SS)  ,C  5, 

With  u"  calc.iid.^«j  “'-X.nT).  Wg  3. 


«‘->»ted,  »urj  .UJ 


“<  ^ u"(x), 


X € ® 


*^<  ^ S(x,  t) , X € 35) 


We  solve  the 


e (0,  T]  . 


(1.3.9) 


equation 


- - 8oo„d.x,  eo„du. 

^ ‘conditions  (1.3.3b) 


(1.3.10b) 


Hoc)  la  given  hy  n 5 2^  t 

‘celc^lleTfr^';^.^*  ;<  ;•  Fina"?iy.  \ 

in  the  .igo- 

■Igorlthm  by  two  bo^n/i  ***  replace  (i  3 J«v’  *"  i"  • method  of 
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u”'*’^(x)  ~ h(x,T) 

(1.3.11) 

where 

h +34  h =0,  x€d  t>0, 

t X 

(1.3.12a) 

h(x,0)  = Ig(x,T,*),  X ^ I!  , 

(1.3. L2b) 

h(x,t)  = u(x,nr+t),  x g SS  , t > 0 , 

(1.3.12c) 

hx(x.t)  “ u^(x,t+nT),  X € t > 0, 

(1.3. 12d) 

and  g 
takes 

satisfies  (1.3.6).  In  this  case  the  solution  of  (1. 
place  through  a sequence  of  convective,  interactive. 

3.7) 

and  dis- 

perslve  processes. 

The  number  u"  Is  generally  dependent  on  n In  (1.3.9),  since 
^ priori  bounds  on  u*^  are  not  obtainable  as  easily  as  for  the 
problems  previously  considered.  In  general  solutions  of  (1.3.7) 
obey  neither  a maximum  nor  a minimum  principle  (References  9 and  12). 
This  is  related  to  the  fact  that,  whereas  each  of  the  components 
(1.3.5)  and  (1.3.6)  of  the  algorithm  for  Burgers'  equation  possesses 
a maximum  and  minimum  principle,  the  same  is  not  true  for  (1.3.12). 

In  the  one-dimensional  case  with  S)  ■ , we  get  maximum  and  minimum 

principles  in  one  direction  only,  in  the  sense  that  if  h(x,0)  = 0 
for  X > Xq , then  h(x,t)  lies  in  the  range  of  h(x,0)  for  x > Xq  and 
t > 0. 


4.  Equations  with  Non-Constant  Coefficients  and  Systems  of  Equations 

Let  ua  restrict  ourselves  to  the  case  where  the  governing 
partial  differential  equation  holds  over  all  space,  so  that  we  can 
focus  attention  on  solution  of  the  initial  value  problem,  and  not 
on  satisfaction  of  the  boundary  conditions.  Taking  u^  ■■  0 in 
(I. I. II)  without  any  essential  loss  of  generality,  we  note  that 
the  algorithm  (1. 1.12)-(I.l. 14)  is  nothing  more  nor  less  than  an 
approximate  solution  of  (1.1.2),  with  initial  condition  (l.l.lb), 
racast  in  the  following  form: 
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u(x,t)  » Urn  u(x,t,0f),  X € R\  t > 0,  (1.4.1a) 

at® 


Uo(x,t)  = 0,  X € , t > 0 , 


(1.4.1b) 


u(x,0,Q')  = min(uo(x),0f),  x € R^ , 0 s a < ® , (1.4.1c) 


(^)  - • 0. 


6 r" 


t>0,  OSjf<oo. 

(1.4. Id) 


In  particular,  (1.4. Id)  is  obtained  by  differentiating  (1.1.2a) 
with  respect  to  or,  and  the  intermediate  solutions  u(of)  are  chosen 
in  such  a way  that  either  the  characteristics 


dx 

dt 


V(u(Qf)) 


(1.A.2) 


of  (1.4. Id)  are  straight  lines,  because  u(a)  and  v(u(a))  are  con- 
stant on  them  (u(o)  - a),  or  else  = 0 on  them.  Thus,  we 

may  view  the  algorithm  (1. 1. 12)-(1. 1. 14)  as  a way  of  building  up 
the  solution  u(x,t),  starting  with  the  trivial  solution  nj(x,t), 
by  means  of  a continuum  of  perturbations,  dependent  on  the  param- 
eter a. 


Similarly,  the  algorithms  of  sections  1.2  and  1.3  may  be 
considered  to  be  obtained  by  formal  differentiation  of  the  equa- 
tions (1.2.1a),  (1.3.1a),  and  (1.3.7a),  respectively,  with  respect 
to  a.  H»e  algorithm  (1.2.12),  (1.3.9),  (1.2.13a),  (1.3.6) , (1.3. 12) , 
and  (1.3.11)  thus  appears  as  an  approximate  solution  to  (1.3.7)  by 
means  of  a continuum  of  perturbations  on  the  solution  u(x,t)  s 0, 
but  it  should  be  emphasized  that  this  is  not  the  same  as  a solu- 
tion built  up  out  of  a formal  expansion  of  the  solution  u(x,t,c) 
of 


uj(x,t,«)  + e V • F(u(x,t,e))  + p A u^(x,t,e)  - 0 


(1.4.3) 
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in  powers  of  e,  a procedure  which  has  been  shown  to  fail  to  yield 
solitary  waves  (Reference  10). 

In  this  vein  we  might  approach  the  problem  of  solving 

G(u,u^,Vu,x, t)  - 0 , (1.4.4) 


subject  to  given  initial  conditions,  as  a problem  of  determining 
a family  of  functions  u(x,t;of),  each  satisfying  (1.4.4),  with 
initial  conditions  dependent  on  or.  By  differentiating  (1.4.4) 
with  respect  to  or,  we  get  the  linear  equation 


‘^u(a) 


du(g) 

9a 


+ G 


Uj.(a) 


Su(a) 


V = 0 . 

oa 

(1.4.5) 


The  characteristics  of  this  equation  are 


dt  ^ ^ 

ds  “u(.(a)’  ds  ”vu(a) 


(1.4.6a) 


and  on  them 

oa 


satisfies 


ds 


+ G - 

9a  ^ ‘'u(a)  9a 


0 . 


(1.4.6b) 


If  X E r",  we  have  n + 2 equations  (1.4.6),  and  the  coef- 
ficients of  these  equations  generally  depend  on  the  2n  + 3 quanti- 
ties X,  t,  u(a),  u^(a),  and  7u(a).  If  we  desire  these  equations 
to  have  constant  coefficients,  for  each  a,  we  can  satisfy  that 
requirement  If  we  pick  a as  a measure  on  the  Borel  sets  In  r!° 
where  m ■ n -f  2,  and  think  of  u(a)  as  a corresponding  signed  mea- 
sure, where  Is  now  to  be  thought  of  as  a Radon-Nikodym 

derivative  (Reference  7).  As  a practical  means  of  solving  (1.4.4), 
this  procedure  docs  not  seem  to  be  generally  promising;  we  indi- 
cate some  more  useful  approaches  below. 
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Higher-order  equations  than  (1.4,4)  can  be  solved.  As  an 
example,  one  might  solve 


Uj.  - uu^^,  u(x,0)  > 0 , 


(1.4.7a) 


by  solving  at  each  time 


t XX 


and  choosing  u(or)  so  that  either  = 0 or  u(a)  = a = constant. 

For  the  equation 


+ XU  =0  , 

t X ’ 


(1.4.8a) 


we  get 


(1.4.8b) 


and  we  want  to  choose  u(a)  so  that  either  ■ 0 or  x ••  of  = constant. 

Note  that  this  equation  is  not  in  conservation  form.  A conservative 
equation  is  obtained  by  writing  u'  > u^  and 


(1.4.9a) 


(IsV 


(1.4.9b) 


For  the  equation 


u^  + xuu  ■ 0 , 

t X 


(1.4.10a) 
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Since  quite  general  nonlinear  systems  of  equations  may  be 
reduced  to  first  order  quasilinear  systems  which  are  homogeneous 
in  the  derivatives  and  do  not  involve  the  independent  variables 
explicitly  (Reference  6),  we  may  think  of  the  equations  above  in 
which  the  independent  variables  appear  explicitly  as  special  cases 
of  homogeneous  quasilinear  first  order  systems  of  equations.  As 
as  example  of  a system  of  the  latter  genre,  consider  the  equations 


G^(u^^\u^®\u^^\u^^^7u^^\vu^®^)  = 0,  i=l,2  , 

(1.4.14) 

.(  = ) 


u<=>.  Vu<^>,  and  Vu' 


where  the  depend  linearly  on  u^ 
the  system  (1.4.14)  is  hyperbolic,  we  may  pick  of  as  a two-dimen- 


If 


sional  measure  such  that  we  either  have 


pici 

uU) 


or  u'  '(a) 


. (3) 

constant  and  u'  (a) 


9of 

constant. 


0 and 


au(_") 


It  should  be  emphasized  that  such  a formal  scheme  to  ob- 
tain solutions  says  nothing  about  the  well-posedness  of  the  origi- 
nal problem,  unless  we  know  something  about  the  well-posedness  of 

all  the  linear  problems  which  are  solved  for  ^ and 


example,  the  system 


For 


(a)  3 / (a)  u^^^) 

“t  ■ ^ « ) - 0 , 


(1.4.15) 


is  only  hyperbolic  if 


Ve  + e / - 4 u'  ''  u'  ' e >0.  (1.4. 


16) 
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CHAPTER  TWO 

convergence  of  the  algorithm  for  a special  case 


The  principal  result  of  this  chapter  will  be  a proof  of 
convergence  of  the  algorithm  (1. 1. 12)-(1. 1. 14)  in  the  case  when 
Xq  ■*  ■“>  +“.  To  carry  out  the  proof,  we  will  have  to  suita- 

bly restrict  the  initial  function  Uq(x)  in  (1.1.1b).  In  the  first 
section  of  this  chapter  we  discuss  suitable  function  spaces  for 
the  study  of  convergence  and  introduce  some  operators  which  ex- 
pedite the  analysis  of  the  algebraic  structure  of  the  algorithm 

(1. 1. 12) -(1. 1. 14) . We  will  be  able  to  reduce  the  study  of  con- 
vergence to  the  case  when  the  initial  function  Uq(x)  is  of  a 
special  type.  In  the  second  section  we  present  an  error  esti- 
mate, and  in  order  to  get  a particularly  sharp  result  we  will 
require  that  the  "velocity"  v(u)  in  (1.1,1a)  be  a monotonic  func- 
tion of  u.  Other  cases  may  be  treated,  and  we  will  show  what 
sorts  of  error  estimates  may  be  expected  for  the  algorithm 

(1. 1. 12) -(1. 1. 14)  in  those  cases. 


1 • Function  Spaces  and  Moaotonlc  Operators 

It  is  clear  that  we  will  need  some  property  of  continuity 
for  the  function  uq(x),  and  cannot  require,  for  example,  that  uq(x) 
just  be  bounded  measurable.  For  instance,  if  we  had 


Uo(x) 


2 2n  To  ^ X < (2n+l)To 

, n=0,+l,±2, • • • , 

0 (2n+l)ToS  x < 2(n+l)To  (2.1.1) 


the  algorithm  (1, 1, 12)-(1. 1. 14)  will  yield  markedly  different  re- 
sults according  to  whether  we  pick  T » To  or  t = To/2. 

We  shall  require  Uq(x)  to  be  almost  uniformly  continuous 
according  to  the  following  definition. 

Definition:  A function  g:  Rl''  -»  is  called  almost  uni- 

formly continuous  (a.u.c.)  if  for  every  a and  b such  that  -<»  < a < 
b < » and  e > 0,  a 6 - 6(c,a,b)  with 


IgCx+Tli  (x)6)  - g(x+Tja(x)5)|dx  < e (2.1.2a) 
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for  all  Tli(x),  11a (x)  such  that 

|t1i(x)1  s 1,  ha(x)l  s 1 . (2.1.2b) 

A more  general  class  of  Uo(x)  will  be  mentioned  later. 

At  this  point,  let  us  recall  the  operators  D and  I intro- 
duced in  (1.1. 17) : 

a < u^ 

u^  s O'  s u(.)  (2.1.3a) 

a > u(  • ) 

> 

IF(.,.)=u<+J  f(-,Qf)do  , (2.1.3b) 

*00 

where 

^ Uo(x),  X € . (2.1.3c) 

Also,  recall  the  definition  of  V in  (1.1.26)  and  (1.1.28). 

Lemma  2.1.1;  Let  u"(x)  be  given  by  (1. 1. 12)-(1. 1. lA) . 

If 

j*  lvjo(x-K)i  (x)6)  - lb  (x-PTla(x)6)ldx  < g (2.1.4) 

a 

for  all  111  and  T)a  satisfying  (2.1.2b),  then 
b-nTV 

/ |u”(x+T)i  (x)6)  - u‘'(x+'n,(x)6)|dx  < G (2.1.5) 

a+mV 

for  all  Such  , f]a . In  other  words,  ug  a.u.c.  =*  u”  a.u.c. 


D(u(*)  ;Qf)  = { 1 
0 


and 
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Proof:  It  is  obviously  sufficient  to  prove  (2,1.5)  for 

the  case  n = 1.  From  (1.1.14)  and  (1.1.15), 


1 u^  (x+\  (x) 6 ) -u^  (x+Tla  (x) 6)  1 s;  J j f°  (x-Hf]i  (x) 6 ,t  ,a) - f®  (x+T]s  (x) 6 ,t  ,a)  [ do 


= j'lf°(x+Tli(x)6-Tv(a),0,o)-f°(x+Tl2(x)6-Tv(o),0,o)ldo  . (2.1.6) 


Then,  by  (1.1.13b)  and  (1.1.12), 


! 

i 


b-VT 

J 


a-fVT 


1 (x-PTli  (x)6)-u‘‘  (x-PTla  (x)6)  1 dx 


® b-VT 

^ J*  I I (x+Tli  (x)6-tv(o)  ,0,o)-f°  (x+Hs  (x)6-tv(o)  ,0,o)  1 dx  do 

-00  a+VT 


00  b-Vr-Tv(o) 

1 I I (x+71i  (x-Kv(o))  ,0,o)-f°  (x+T^  (x+Tv(o))  ,0,o)  jdx  do 

-00  a4VT-TV(o) 


b 00 

^ J J*  lf°(x+TH  (x+Tv(o)),0,o)-f°(x+Tls(x+Tv(o)),0,o)ldo  dx 

a -oD 


b 

sj  / sup  uo(y)  - inf  Uo(y)'\dx 
a 'y€[x-6,x+6]  y€[x-6,x+6]  / 


b 

|uo(x+Tl3(x)6)-uo(x+T)4(x)6)ldx  + 2 6(b-a)  , (2.1.7) 


where  we  pick  Tb  (x)  and  TU(x)  to  satisfy  the  same  requirement  (2.1.2b) 
as  TH  and  T)a  ■ and 


uo(x41t(x)6)  a sup  Uo(y)  - 0 , 
y€[x-6  ,x+6] 

Uo(x41U(x)6)  S Uo(y)  + e , 

y^Lx-o  ,x+6] 


(2.1.8) 
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with  0 > 0.  Letting  6 I 0 and  using  (2.1.4),  we  obtain  (2.1.5) 
for  n - 1. 

Before  proceeding  further,  let  us  introduce  operators  E", 
where  T]  is  a function  of  a:  ^ 


f)  (X,*)  = sup  f(y,*) 

' y€i(x.Ti) 


where 


(2^  f)  (X,*)  = inf  f(y,-)  , (2.1.9b) 

y^i(x,'n) 


I(x,Tl)  = [uin(x,x+Tl(.)),  max(x,x+7K-))l  . (2.1.9c) 


The^following  properties  can  readily  be  verified  for  the 
operators  E“  : 


E^  E^  - E^+  E'*’-  . 

Til  ^8  \ +Tla  -Ffla  ’ 


where 


E"  E"  - E’_,  E". 

y\a  Til +713 


71  (.)  s max(71(-),0) 


71  (•)  = min(71(*)i0) 


(2.1.11) 


In  particular,  it  follows  from  (2.1.10)  that 


EiJ;  E^  - E'^  E"*" 
TTi  us  7V»  7H 


(2.1.12a) 
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and  if  TIiTla  i:  0, 


E"  » z"  e" 

Til  Tla  % 


\ Tb 


lb  ' %-ni2 


Also 


and 


If  |ti1  s 1 6|  , we  get 


and 


E^  ZL  £ E"  E"*" 
Til  MS  Tls  Til 


^’■n  K ^ 1 
-T)  7) 


T1  -T1 


E"  e;!'  e:  - e: 

-T1  T1  6 6 


E’^  Et  - Et 
n -T1  6 6 


In  terras  of  the  E's  we  can  define 


“ E^  E^ 

11  T1  ^-Tl 


(2.1.12b) 


(2.1.13a) 


(2.1.13b) 


(2.1.14) 


(2.1.15b) 

(2.1.15b) 


(2.1.16a) 


(2.1.16b) 


(2.1.17) 


\ 
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Given  a function  g(x)  and  a number  6 > 0 independent  of  a,  we  can 
also  define 


(Sg  g)  (x)  - (ISg  Dg)  (x)  , (2.1.18) 


In  other  words 


(Sg  g)  (x)  - sup  g(y)  , (2.1.19a) 

y€[x-6  ,x+6] 


(Sg  g)  (x)  » inf  g(y)  . (2.1.19b) 

y€[x-6,x+6] 


Finally,  let  us  introduce  the  shift  operator 


(T^f)(x,*)  » f(x-!l(.),')  . (2.1.20) 


The  following  relations  among  these  operators  can  easily  be  estab- 
lished : 


(2.1.21a) 


%/2  - S, 


-T)  -T1  Tl/2  * 


(2.1.21b) 


and 


=V2 


’ll  ®T|/2 

We  get  conmutativity  between  the  E's  and  T: 


(2.1.21c) 


(2.1.22) 
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and  we  note 


± + 
i:  . - T-  I ■ 

-T\  T1  T1 


(2.1.23) 


so  that 


” ’■-n  ^-all 


(2.1.24) 


and 


All  the  operators  Introduced  are  monotonic.  That  is, 


8i  ^ ga  =*  Ogi  i Oga , 0 * D,  S.,  or  s"  (2.1.25a) 


fi  ^ fa  =»  Ofi  i Ofa.  0 - I,  T^,  sijj,  s’. 


From  the  definitions  it  follows  that 


> s;  I 


DsJ  - D 


or  Z. 


'1^ 


(2.1.25b) 


(2.1.26a) 

(2.1.26b) 

(2.1.26c) 


and 


ISg  s Sj  I 


(2.1.27a) 


DS.  '■  S.  D 

0 u 


(2.1.27b) 


(2.1.27c) 
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Also, 


and 


(2.1.28a) 


S.  D T D s S,  D 
6 T]  6 


(2.1.28b) 


where  ItKot)!  s 6 for  all  a. 


In  terms  of  these  operators,  with  the  help  of  (1.1.15),  we 
can  state  the  algorithm  (1. 1. 12)-(1. 1. 14)  as 


n+i  n 

u - IT _ Du 
VT 


(2.1.29) 


Thus,  the  algorithm  proceeds  entirely  through  the  application  of 
monotonic  operators.  If 

Uo(x)  ^ So(x)  , (2.1.30a) 


and  we  denote  the  respective  approximate  solutions  by  u"(x)  and 
u*'(x)  (with  the  same  u^  for  Uq(x)  and  Uo(x))  we  get 


u"(x)  2 5"(x)  . (2.1.30b) 


Thia  sort  of  monotonicity  result  also  holds  for  the  problem  (1.1.5) 
posed  on  a finite  domain,  provided  we  have  the  appropriate  monoto- 
niclty  required  of  the  corresponding  boundary  data.  The  result 
also  holds  for  the  algorithm  (1.2. 12)-(1.2. 14)  used  to  solve  the 
problera  (1.2.1),  and  for  the  algorithm  used  in  the  solution  of 

(1.3.1) .  Assuming  convergence  of  the  algorithms  for  these  prob- 
lems, we  conclude  that  the  monotonicity  carries  over  to  the  exact 
solutions  c'  (1.1,5),  (1.2.1),  and  (1.3.1)  generated  by  the  algo- 
rithois.  At  this  juncture  we  should  reemphasize  the  point  made  in 
•®ctlon  1.1,  that  the  solution  to  (1.1.1)  computed  by  our  proce- 
<lura  is  not  necessarily  the  desired  one  In  terms  of  any  underlying 
physical  rationale.  For  example,  the  one-dimensional  flot^  of  an 
invlscld  prassureless  fluid  is  governed  by  an  equation  of  type 

(1.1.1)  with  u the  velocity  and  v(u)  ■ u.  However,  the  underly- 
ing physical  laws  are  conservation  of  mass  and  momentum,  and  u is 
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really  coupled  to  a density  field  p.  For  general  initial  density 
fields,  we  do  not  have  monotonic  dependence  of  u(x,t)  on  u(x,0). 
Hence  any  algorithm  with  that  property  will  have  to  be  in  error, 
in  general,  in  a later  part  of  this  work  we  will  present  a phy- 
sically correct  algorithm  to  solve 


u.  + uu  “0 

t X 


when  u is  the  velocity  of  a one-dimensional  flow. 

Returning  to  the  operators  Sj  , we  see  that 


^6  Sj  ^ 1 , 
^6  ° ^6 


(2.1.31) 


The  operator  Sj  "flattens  out"  all  peaks  of  a function  g(x)  by 
removing  portions  of  width  < 26.  Similarly,  s"  "flattens  out" 
the  troughs  by  filling  in  the  portions  of  width  <'*26. 

Upon  application  of  equations  (2.1.26),  (2.1.27),  and 
(2.1.31),  we  see  that  the  following  inequalities  can  be  estab- 
lished by  induction  in  n: 


sJa^D)"  uo  ^ (SjsJjS-lT^D)"  s;  sj  Uo 


2 S 


6 si  Uo 


(2.1.32) 


significance  of  (2.1.32)  is  that  it  gives  bounds  on  the  quan- 
tities w"  computed  by  the  algorithm 


■ ®6  ^6  ®6 


(2.1.33a) 
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- s"  s’*",  s"  IT  D w"  . (2.1.33b) 

6 a 6 6 VT 

When  Uq  is  a.u.c.,  we  can  use  these  bounds  to  show  that  w*^  is 
close  to  u*^,  as  we  shall  do  below.  Assuming  that  this  is  the 
case,  we  see  then  that  we  may  equally  well  use  the  algorithm 
(2.1.33),  as  use  the  original  algorithm  (1. 1. 12)-( 1. 1. 14)  (or 
(2.1.29))  to  solve  the  problem.  Furthermore,  each  of  the  quan- 
tities w",  n 2 0,  has  the  property  that  all  peaks  and  troughs  are 
of  width  2 26.  This  fact  will  permit  us,  in  the  convergence  proof 
given  in  the  next  section,  to  restrict  our  attention  to  portions 
of  the  profile  w (x)  which  are  monotonically  non-increasing  or  non- 
decreasing. That  is,  if  for  example  (^o»?i)  is  the  interval  be- 
tween a flattened  peak  and  trough  of  Uq , then  for  0 ^ nr  ^ 26/V 
the  function  w”(x),  x€  [§o>?ili  will  be  either  non-increasing 
or  non-decreasing,  according  as  v/’  is.  It  is  clear  that  if  we 
can  establish  convergence  and  error  estimates  as  t 1 0 for  such 
intervals  (x,t)  € X [0,  26/v],  then  we  can  do  so  for  all 

(x,t)  € (-”,<*>)  X [0,«). 

To  show  that  w"  is  close  to  u"  if  6 is  sufficiently  small, 
we  use  (2.1.32)  in  the  form 


^ J.  ^ 

^ W ^ 


(2.1.34) 


where 


6P  -S^uo 


(2.1.35a) 


and 


0n+i  __  _ on 

- IT  D 9 , 

VT  ’ 


“ S7  uo 


xO+i  __  ..  j n 

♦ - IT  , D t 

VT  ^ 


(2.1.35b) 

(2.1.36a) 

(2.1.36b) 


Because  of  the  monotonicity  result  (2.1.30),  we  see  from  (2.1.35) 
and  (2.1.36)  that 


e"  2 


j-  ‘ 

U 2 t 


(2.1.37) 
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and  hence 


lw"-u"l  s e"  - . (2.1.38) 


Because  of  the  monotonlc  dependence  of  6*^  and  t|("  on  the 
initial  data  9^  and  , and  the  fact,  which  follows  from  the  de- 
finition of  V,  that  all  non-zero  contributions  to  6"  and  i|r"  in 
[a+nTV,b-nTV]  have  originated  as  contributions  to  6^  and  \|(°  some 
where  in  [a,b],  we  get 


b-nTV 

; 


a+ntV 


b 

(0%")dx  s r (0P-^O)dx 
a 


(2.1.39) 


From  (2.1.38),  (2.1.39),  and  the  initial  conditions  (2.1.35a)  and 
(2.1.36a),  we  get  the  following  result. 

^ Leonta  2.1.2;  If  Uq(x)  is  a.u.c.  and  satisfies  (2.1.2)  and 
if  w is  given  by  the  algorithm  (2.1.33)  and  u"  is  given  by  (1.1.12)- 
(1.1.14),  then 


b-nTV 

I |w"(x)-u"(x)ldx  s e(6)  , (2.1.40) 

a+nrV 


where  we  denote  the  right  hand  side  of  (2.1.2a)  by  e(6,a,b). 

With  the  aid  of  the  operators  introduced,  we  can  also  find 
the  effect  on  the  approximate  solution  u"(x)  of  (1.1.1)  due  to  a 
perturbation  of  the  "velocity"  v(u),  when  Uq(x)  Is  a.u.c.  On  the 
assunption  that  the  approximate  solution  converges  to  an  exact  so- 
lution, the  same  error  estimate  will  carry  over  to  the  exact  solu- 
tion (with  nr  replaced  by  t and  u"(x)  replaced  by  u(x,t)). 

Lemma  2.1.3:  In  the  case  Xq  "*  ■ ®i  X;  "•  + ®,  let  u”(x) 

be  generated  from  uq(x)  by  equations  (1. 1. 12) -(1. 1. 14) , and  let 
u be  generated  from  Uq(x)  by  the  same  equations  with  v(o)  re- 
placed by  v(of).  Suppose  Uq(x)  is  a.u.c.  and  satisfies  (2.1.2). 

Use  the  notation  e(6)  In  (2.1.40).  Suppose  lv(a)l  s V and 
lv(a)l  < V for  all  er  € [u  ,U],  where  U is  given  In  (1.1.26).  Then, 
with  |v(er)-v(or)l  < v'  for  all  a 6 [u  ,U],  we  have 
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b-ntV 

J*  |u"(x)-u"(x)ldx  s e(nTV')  . (2.1.41) 

a+nTV 


I 

I 

I 

I 

I 

I 

I 

I 


I 

I 

I 


Proof:  Using  (2.1.29),  we  can  write 


u"  - (IT^^D)"  uo 


(2.1.42a) 


S"  - (IT^D)"  Uo 


(2.1.42b) 


Now,  we  can  establish  inductively  that 

Cv'  . (2.1.43) 

The  case  n - 0 is  obvious.  To  get  from  n to  n + 1,  we  note  that 


^ ISJ/,  V D u"  ^ IT^  D (IT^D)"  Uo 

upon  repeated  application  of  (2.1.26).  Similarly,  we  can  show  that 


~n 


VT 


Let 


®°n  “ 


VT  n* 


®nTV' 

(2.1.44) 

n ^nxV'  ^ ’ 

(2.1.45a) 

" - (IT^D)"  • 

n VT  ' ^n 

(2.1.45b) 

That  is,  0”  and  ara  the  functions  generated  at  the  nth  time 
step  for  tne  unperturbed  velocity  v(a),  from  the  initial  functions 
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0^  and  , respectively.  By  reasoning  similar  to  that  which  led 
to  (2.1.39), 


b-ntV  b 

a+nxv  a 

The  monotonic  dependence  on  the  initial  data  also  implies 


. <2.1.47) 


From  (2.1.47),  (2.1.43),  (2.1.44),  (2.1.45),  and  (2.1.46),  we  get 


b-nTV  b 

' ll”(x)-u"(x)ldx  s J"  (S*  Uo-S"  . u„)dx  , 

“ (2.1.48) 


and  this  leads  directly  to  (2.1.41),  upon  use  of  the  almost  uniform 
continuity  of  uq  and  (2.1.2). 

Thus  if  v(or)  is  piecewise  continuous,  we  can  replace  it  by 
a piecewise  constant  function  v with  Iv-vj  as  small  as  we  desire, 
and  for  a.u.c.  initial  data  get  as  small  a change  in  the  solution 
un(x)  as  we  desire  (in  the  sense). 

By  an  argument  similar  to  that  used  to  establish  lemna 
2.1.3,  we  can  prove  the  continuous  dependence,  in  the  sense, 
of  the  functions  generated  by  (1. 1. 12)-( 1. 1. 14)  on  the  "time". 

We  again  use  monotonicity  and  note  that 


There  follow 


m 


& u 


xf 

mrV 


(2.1.49) 


b-nrV 

J |u  (x)-u”'*®(x)ldx  s e(mrV)  . (2.1.50) 

a+nTV 
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A somewhat  broader  notion  of  continuity  is  that  of  "almost 
continuity"  (a.c.)*  of  a function  g,  by  which  we  mean  that,  for 
any  e > 0,  we  can  find  a set  Sg  of  Lebesgue  measure  < e such 
that  g is  uniformly  continuous  on  [a,b]  - Sg.  We  will  find  this 
type  of  continuity,  suitably  generalized  to  flows  with  other  meas- 
ures, to  be  especially  Important  in  the  treatment  of  actual  physi- 
cal flows.  Accordingly,  there  will  be  a more  complete  discussion 
of  almost  continuity  in  a later  report  in  this  work.  Essentially, 
instead  of  removing  a set  of  Lebesgue  measure  < e from  the  domain, 
we  will  want  to  remove  sets  of  small  mass  or  action,  depending  on 
the  precise  type  of  continuity  considered,  from  various  domains  of 
independent  variables.  Further  on,  when  we  treat  the  flows  in  a 
stochastic  framework,  still  further  enlargement  of  the  notion  of 
almost  continuity  will  be  in  order.  For  the  present,  let  us 
merely  note  that  if  g is  a.c.,  then  g is  measurable,  but  the 
reverse  implication  does  not  hold.  In  addition,  it  follows  from 
the  definitions  that  if  g is  a.u.c.,  then  it  is  a.c.  Ihe  converse 
is  not  necessarily  true,  unless  the  set  Sg  is  dense  on  a set  in 
[a,b]  whose  measure  1 0 as  e i 0.  For  example,  the  converse  is 
not  true  if  g is  1 at  all  rational  numbers  in  [a,b]  and  0 else- 
where. Then  g is  a.c.  but  not  a.u.c.  With  a little  reflection, 
we  see  that  an  analogue  of  lemma  2.1.1  can  be  proven  when  Uq  is 
a.c.:  Uo&«c.  ^ u^a.c.  Similarly,  the  continuous  dependence  of 

the  functions  u^(x)  generated  by  the  algorithm  (1. 1. 12)-(1. 1. 14) 
on  the  "velocity"  v(a) , stated  by  lemma  2.1.3  for  a.u.c.  initial 
data  uq(x),  can  be  proven  for  a.c.  initial  data,  and  the  suita- 
bility of  replacing  the  algorithm  (1. 1. 12)-(1. 1. 14)  by  the  algo- 
rithm (2.1.33),  stated  in  leoina  2.1.2  for  a.u.c.  initial  data, 
holds  for  a.c.  initial  data. 

The  types  of  continuity  we  have  been  considering  in  this 
section  are  weaker  than  that  required  in  section  1.1  to  show  that 
the  right  hand  sides  of  (1.1.33)  and  (1.1.37),  for  0 < t & t,  are 
0(1)  and  o(t),  respectively,  as  t 1 0.  It  would  be  tempting  to 
conclude  that  Uo(x)  a.u.c.  or  a.c.  implies  that  u(xo,C)  is  a.u.c. 
or  a.c.  in  time  for  each  Xq.  This  is  not  the  case,  however,  as 
we  can  see  from  the  following  example.  In  the  limit  i 0,  and 
with  v(or)  ■ a,  choose  the  initial  profile  near  x - 0 to  satisfy 
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Uo(x) 


U - 0x 

0 > x > -UY 

-U 

0 < X < UY 

U(l-I^Y) 

-UY  > x > -2UY 

-U  - 3(x-UY) 

UY  < X < 2Uy 

(2.1.51) 


I^'itially  there  is  a jump  2U  at  x = 0,  and  it  is  easy  to  see  that, 
if  s(t)  is  the  location  at  time  t of  the  discontinuity  in  u(x,t), 

^ (0)  = 0.  After  a time  2y,  we  get  s(2y)  = U,  u(s(2y)',2y)  = 
U(l+eY),  u(8(2Y)t2Y)  = -U(l-lfY),  and  ^ (2y)  = 0.  We  can  now  ad- 
just Uq(x)  so  that  s(t)  moves  back  when  2y  < t < 4y,  forward  again 
when  4y  < t < 6y,  |nd  so  on.  After  a time  t,  in  the  limit  ^ i 0, 
we  get  a jump  2Ue^  at  a position  s(t)  which  oscillates  between 


if  )u  ef  u . 

Points  in  the  interval  [0,  U]  are  Included  in  all  the  oscilla- 

tions. In  particular,  we  cannot  make  any  statement  about  the  a.c. 
or  a.u.c.  of  u(0,t)  which  will  be  valid  Independent  of  the  choice 
of  Y.  although  the  initial  profile  Uo(x)  will  be  a.u.c.,  accord- 
ing to  the  definition  (2.1.2),  with  e(6)  independent  of  y.  The 
values  of  x for  which  we  get  uniformly  large  variations  in  u(x,t) 
when  t changes  by  an  amount  Y independent  of  x will,  of  course, 
have  Lebesgue  measure  which  goes  to  0 as  Y ^ 0,  by  the  result 
(2.1.50). 

(The  fact  that  we  cannot  get  a priori  that  u(xo,t)  is 
a.u.c.  or  a.c.  for  each  Xq  is  related  to  our  ability  to  control 
the  free  boundary  problem.  In  particular,  for  the  case  v(a)  - a, 
if  we  are  given  quantities  J(t),  s(t)  for  0 s t « T with  the  follow- 
ing properties: 


J(0)  2 0 , 


(2.1.52a) 


^ (t)  bounded,  non-negative,  0 s t s T , 


(2.1.52b) 
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I d^8 


(t)l 


2 dt 


0 s t S T 


(2.1.52c) 


then  there  is  an  initial  profile  Uo(x),  which  is  uniqut  id  differ- 
entiable in  the  Intervals 

s(T)  - T J(T)  + ^ (T)J  ix<0  (2.1.53a) 

and 

0 < X 8(T)  + J(T)  - (T)j  , (2.1.53b) 

and  which  satisfies 

X + Uo(x)T  < s(T)  for  x < s(T)  - t[j  J(T)  + — (T)j  , 

(2.1.54a) 

X + Uo(x)T  > s(T)  for  x > s(T)  + t|j  J(T)  - ^ (T)j  , 

(2.1.54b) 


such  that  the  solution  of  (1.1.2),  (1.1.1b)  has  the  following  prop 
erty : 


u(8  (t),t)  - u(s''’(t),t)  - J(t),  0 ?;  t s T . (2.1.55)) 


It  follows  from  the  preceding  discussion  that  we  should 
only  expect  equations  (1.1.6)  and  (1.1.7)  to  be  satisfied  by  the 
function  u(x,t)  “ {(u”(x))^_^^^} , when  it  exists,  in  some 

integrated  sense  (in  time). 

If  we  average  equation  (1.1,23)  over  an  interval  of  values 
of  a and  b,  so  that  the  left  hand  side  becomes 


h b-m 

■Jh*  J*  X f‘'(x,t,o)dx  dl^  , 
-h  a-H] 
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then  we  will  obtain  averaged  versions  of  (1,1.6)  and  (1.1.7)  which 
will  hold  in  the  limit  t I 0,  by  two  applications  of  the  continuity 
result  (2.1.50),  with  a + nTV,  b - nTV  replaced  by  a-h,  a+h  and 
b-h,  b+h,  respectively. 

Although  we  have  not  discussed  the  case  of  the  boundary 
value  problem  on  the  finite  interval  [xo,x^]  in  this  section,  we 
may  expect,  from  the  discussion  following  (1.1.33),  that  the  ap- 
pi^opriate  data  to  specify  at  the  boundary  points  Xq  and  Xj  are 
u (Xo,t)  and  u”(x^,t)  a.u.c.  or  a.c.  in  L. 


2 . Convergence  and  Error  Estimate 

The  principal  result  of  this  section  is  that,  when  the  veloc- 
ity v(u)  is  a monotonically  non-increasing  or  non-decreasing  func- 
tion of  u,  then  the  thickness  of  "shocks"  generated  by  the  algo- 
rithm ( 1 . 1 . 12) -( 1 . 1 . 14)  is  never  more  than  2Vt , where  V is  given 
by  (1.1.26)  and  (1.1.28).  Furthermore,  the  shocks  so  generated 
are  within  the  distance  2vt  of  their  position  in  the  exact  solu- 
tion. In  order,  however,  to  indicate  how  one  may  obtain  other 
types  of  error  estimates  which  might  be  more  useful  in  certain 
situations,  and  to  show  why  we  will  want  to  restrict  outselves  to 
r/'fotone  V if  we  are  to  get  particularly  sharp  error  bounds,  we 
shill  proceed  initially  in  a somewhat  more  general  vein.  Our  de- 
liberations  will  be  expedited  by  the  results  we  obtained  in  the 
I last  section, 

i 

- So  far  we  have  not  stated  what  assumptions  we  make  on  the 

function  v(*)»  But  it  is  clear  from  (2.1.29)  that,  if  the  operator 
i "I"  in  the  algorithm  is  to  be  meaningful,  then  we  should  at  least 

require  v to  be  measurable.  When  v is  piecewise  continuous  and 
* Uq(x)  is  a.u.c.,  we  h.ive  already  seen,  in  lenina  2.1.3  and  the  dis- 

cussion following  it,  that  we  may  replace  v by  a piecewise  constant 
; function  v.  The  study  of  convergence  of  the  algorithm  (1.1.12)- 

(1.1.14)  for  piecewise  constant  v should  be  a relatively  simple 
natter,  since  we  may  imagine  that  we  are  dealing  with  the  propaga- 
. tion  of  a finite  number  of  Individual  "streams",  each  moving  with  a 

constant  velocity,  with  a specified  law  of  interaction,  and  we  know 
from  Chapter  One  that  the  algorithm  gives  the  exact  solution  when 
arch  of  these  streams  propagates  Independently. 

In  fact,  when  Uq(x)  is  a.u.c.  it  should  be  possible  to 
prove  convergence  when  all  we  require  is  that  v be  bounded  measura- 
ble. For,  by  the  monotonicity  properties  shown  in  the  last  section. 
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we  may  bound  u above  and  below  by  6 and  respectively,  where 


^(x)  = sup  Uo(y),  x6  [m6,(nri-l)6)  , (2.2.1a) 

ye[m6,(m+i)6) 


Sn+i  „ ~n 

0 = IT  D e , 

VT  ’ 


(2.2.1b) 


/"T 

■V 


and 


t°(x)  - Inf  Uo(y),  x€  [m6 , (m+l)6)  , (2.2.2a) 

y€[m6 , (tn+i)6) 


V = IT  D * . 

VT  ' 


(2.2.2b) 


If*  (2^.1)  and  (2.2.2),  m is  an  integer,  - oo  < n < a>.  Clearly 
^ ^ ^ Sg  Mo  and  uo  ^ Sg  , Accordingly,  the  argument 

used  to  prove  lemma  2.1.2  can  be  used,  and  we  find 


b-nTV 

1 6*'(x)  - u”(x)  1 dx  s 

I 

a+niV 

e(6)  , 

b-nTV 

J 

a+nTV 

It  (x)-u"(x)ldx  s 

e(6)  . 

Note  that  the  initial  data  9°  and  in  (2.2.1a)  and  (2.2.2a)  are 
"histograms".  Hence  when  v is  bounded  measurable  we  need  only 
prove  convergence  for  such  initial  profiles.  The  exact  solutionr, 
for  such  profiles  can  be  given  for  such  v.  We  know  that  a "Jump' 
at  X ■ s(t)  from  uq  to  u^  moves  with  velocity 


ds  F(ui)-F(‘k)) 

dt  ” u^ -Uq  ’ (2*2.4) 


where  F is  given  by  (1.1.2b).  In  fact,  when  the  initial  profile 
is  such  a histogram  involving  only  a finite  number  of  values  of  u, 
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say,  {u£,  it  is  clear  that  the  exact  solution  is  un- 

changed by  changes  in  the  velocity  field  v which  leave  differences 
of  values  of  F at  these  values  {uj^}  invariant.  What  sort  of  con- 
vergence result  we  may  expect  for  our  algorithm  for  v bounded  mea- 
surable when  the  initial  profile  is  a jump,  will  be  Indicated  in 
the  sequel. 

By  lemma  2.1.2,  we  may  deal  with  functions  w"  whose  peaks 
and  troughs  are  all  of  width  ^ 26.  As  we  pointed  out  in  the  dis- 
cussion accompanying  the  lemma,  this  permits  us  to  restrict  our- 
selves to  a portion  of  a profile  v^(x)  which  is  monotonically  non- 
decreasing or  non- Increasing.  Without  loss  of  generality,  we  take 
the  latter  case.  Given  such  a portion  of  a profile,  say,  for 
s X S §1 , we  may  let 


Then 


U “ inf  fcf’(x). 


U*"  = sup  y>P  (x) 

>'^[§0.?l] 


(2.2.5) 


lim  v^(x)  ■ U^,  lira  w®  (x)  = U , (2.2.6) 

X i ?0  X t 


and  we  may  imagine  the  initial  profile  v/*  (x)  extended  so  that 


^ (X)  - u‘*‘, 

X < . 

(2.2.7a) 

wP(x)  - U", 

X > §1  . 

(2.2.7b) 

With  this  profile,  we  find  that 

w "IT 

VT 

D w'^'  , 

(2.2.8) 

so  that  effectively  we  are  considering  again  the  functions  gene- 
rated by  the  original  algorithm  ( 1. 1. 12)-(1. 1. 14) , or  (2.1.29). 
Thus,  we  shall  now  write  \P , u"  for  \P  , w". 
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I Before  we  proceed  with  a discussion  of  convergence,  let  us 

t write  down  the  exact  solution  of  (1.1.2)  with  initial  data  u(x,0)  » 

i/’Cx).  The  proof  of  convergence  and  error  estimate  will  then  be  ob- 
tained by  comparing  the  approximate  solutions  u”(x)  with  the  exact 
I solution  u(x,nr). 

I 

; With  a monotonlcally  non-increasing  profile  g(x),  lim  g(x)*U^, 

y-®  g(x)  = U , we  associate  a function  X(g;of)  with  range  the  ex- 

^ - + 

tended  real  line  for  a € [U  ,U  ],  as  follows: 


(1  X i X(g;a) 

(Dg)(x,a)><  . (2.2.9) 

(0  X > X(g;a) 


For  the  solution  of  (1.1.2)  with  u(x,0)  = u°(x),  we  use  the  notation 
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^ Xo(a')dor'  F(ao)-F(a) 

inf  + t 

Q^[u',ao]  Qfo  - a oto  - a 


In  this  way  we  proceed  until  we  have  located  all  the  shocks. 
More  precisely,  given  c > 0 and  P € (U  ,U^],  let 


. jy  Xo(Qf  )do  f(P)-F(Y)  \ 

ala€[U  ,e),  or  - inf  ( + -—r — tl<e  , 

v€[U-,e)^  ^ ^ &-Y  / 

(2.2.13a) 


A (P,t)  * closure  of  A (0,t) 
€ € 


(2.2.13b) 


Define 


“o(0.t)  = inf  [of|o€Ag(3,t)  > 0} 


(2.2.14) 


Of  (e.t)  = Oo(a  (e,t),t)  . 


(2.2.15b) 


Then  the  exact  solution  of  (1.1,2)  with  u(x,0)  - u°(x)  is  given, 
for  3 € (UMT*-],  by 


is  non- Increasing  in  or.  Let  the  infemun  in  (2.2.12)  be  achieved 
for  oq*  We  then  proceed  to  find  the  next  shock  in  the  continuum 
of  shocks.  (Those  who  object  to  the  "rigor"  of  the  argument  may 
prefer  that  we  start  with  an  initial  u(x,0)  given  as  a non-increas- 
ing histogram.  Tlien  the  shocks  will  be  denumerable.)  The  "next" 
shock  will  be  located  at 


It  is  clear  that  cro(3,t)  is  non- increasing  in  p.  Further,  set 


S sup  {ol3€[Qfo(Of,t),Qf]) 

a€[u‘,U+] 


(2.2.15a) 
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x^(3)  - lim 
ei  0 


Xo(a')do' 


O'  (B,t)-a  (3,t)  +€ 


F(Qf+(3.t)-F(a-(B,t)-e) 

— t 

Of  (3,t)-0f  (3,t)+  e 


(2.2.16) 


Returning  to  the  approximate  solutions,  we  see  that,  if 
g(x)  = h(x+6),  then 


(IT^^Dg)(x)  - (IT^Dh)(x+6)  . (2.2.17) 


Accordingly i if  g is  monotonlcally  non-increasing  and  6 i 0,  so 
that  g(x)  i h(x),  we  get  from  the  monotonicity  of  the  operators 
I,  D that 


(IT^Dg)(x)  s (IT^^Dh)(x)  - (I'^^^Dg)(x-6)  , (2.2.18) 


and  Jhus  the  functions  u"(x)  will  be  monotonlcally  non-increasing 
If  u is.  Given  u"(x),  using  the  notation  (2.2.9),  we  define 


X*(of)  - X(u"(.);o)  . (2.2.19) 


In  addition,  we  define  x*(of) , a € [u”,U^],  by 


Note  that 


(T^^Du")(x,a)  - 


1 X s X*(Of) 

0 X > x*(of) 


(2.2.20) 


X*(0f)  + Tv(of) 


(2.2.21) 
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and 


Xj(of)  « xJ(cTi  ) if  a a cYi 


(2.2.22) 


IL  might  seem  that,  if  for  some  n and  a > we  have 
x*(or)  > x*(«i)  (thereby  implying  by  (2.2.21)  and  (2.2.22)  that 
VLO;  > v(,ai)),  we  would  have  1 x^(a)-xS(Oi  ) 1 ^ 2Vt  for  m ^ n. 

Such  a result  could  then  be  used  to  bound  the  width  of  a numeri- 
cal "shock",  or  rapid  change  in  u including  [otj^  ,01],  and  lead  to 
a convergence  result  and  error  estimate.  Our  expectation  is 
wrong,  however,  as  we  may  see  by  considering  the  example 


u°(x) 


1 X < 0 

0 X > 0 


(2.2.23a) 


with  the  velocity  profile 


v(a)  - 


0 

1 - c 

1 


1 < o < 1 


6 > 0 


(2.2.23b) 


The  exact  solution  of  (1.1.2)  with  initial  condition  (2.2.23a)  and 
velocity  profile  (2.2.23b)  is 


u(x,t) 


X < - j e t 
-je  t<x<0 

X > 0 


(2.2.24) 


We  have  not  shown  that  the  approximation  to  (2.2.24)  does  not  have 
numerical  shocks  whose  width  is  bounded  by  a multiple  of  t,  but 
only  chat  a conjecture  which  would  readily  lead  to  a boupd  of  that 
type  is  invalid.  The  violation  of  Che  conjecture  seems  to  have 
something  Co  do  with  Che  non -monotonicity  of  v(a)  in  (2.2.23b). 


1 

i 


I 

I 


I 

I 

I 

1 

I 


I 

I 

1 

1 

] 


1 

3 

i) 
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A more  graphic  demonstration  that  we  cannot  hope  to  have 
numerical  shocks  of  thickness  0(t)  as  t i 0 when  v Is  merely  re- 
quired to  be  measurable  is  given  by  the  following  example.  Sup- 
pose that  v(of)  is  measurable  and  takes  on  the  values  0 and  1, 
such  that  the  measure  of  the  subset  of  any  interval  in  which  v is 
1 is  one  half  the  Lebesgue  measure  of  the  interval,  and  the  measure 
of  the  subset  in  which  v is  0 is  also  one  half  the  Lebesgue  measure 
of  the  interval.  Consider  the  initial  profile  given  in  (2.2.23a). 
We  know  that  the  exact  solution  to  the  problem  is 


u(x,t) 


On  the  other  hand,  we  compute 


X < 0 

0 < X < T , 
X > T 

X < 0 
0 < X < T 

T < X < 2t 
X > 2t 

X < 0 
0 < X < T 

T < X < 2t  , 

2t  < X < 3t 
X > 3t 


(2.2.25) 


(2.2.26a) 


(2.2.26b) 


(2.2.26c) 
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u*(x)  - 


X < 0 
0 < X < T 
T < X < 2t 
2t  < X < 3t 

3t  < X < 4t 
X > 4t 


(2.2.26d) 


and  so  on.  It  is  apparent  that  we  get  a binomial  distribution  for 
the  "front",  which  becomes  a Gaussian  distribution  in  u"(x)  as  nt“. 
Accordingly,  the  width  of  the  front  is  0 (v' nVT ) » 0(VtTV)  where 
t ■ nT . In  this  case  the  time  step  t gives  an  effective  "viscosity" 
0(V®T). 


A natur«"V  conjecture  to  make  is  that  when  v(*)  is  bounded 
measurable  the  width  of  shocks  is  0(^/ trV) . We  do  nothing  further 
in  this  report  to  establish  the  conjecture,  but  it  would  appear  to 
be  susceptible  to  proof  by  probabilistic  methods. 

For  the  remainder  of  this  section  we  focus  on  the  case  where 
V is  monotonic.  Since  we  have  already  required  the  initial  profile 
(x)  to  be  non-increasing  and  we  know  from  section  1.1  that  the 
algorithm  ( 1. 1. 12)-(1. 1. lA)  gives  the  exact  solution  to  the  problem 
in  this  case  when  v la  axmotonically  non-increasing,  we  shall 
consider  the  other  case,  where  v is  monotonlcally  non-decreasing. 
That  is,  we  require 


v(cy)  2 v(aj)  if  O'  2 Qf^  . (2.2.27) 


With  the  assumption  (2.2.27)  and  the  various  simplifica- 
tions of  the  general  initial  value  problem  which  have  been  made, 
we  could  plunge  directly  into  a proof  of  our  central  result.  We 
will  do  this  shortly,  but  first  we  will  consider  some  simpler  re- 
sults which  illustrate  the  ideas  we  use  in  carrying  out  the  proof. 
There  will  be  some  further  discussion  of  the  types  of  situations 
which  BMy  arise  in  implementing  the  algorithm  (1. 1. 12)-(1. 1. 14) , 
in  order  to  better  motivate  the  language  in  which  our  central 
result  is  cast. 
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We  begin  by  verifying,  for  the  case  (2.2.27),  the  conjec- 
ture tentatively  put  forth  after  equation  (2.2.22),  and  shown  to 
be  Invalid  for  general  velocity  profiles  v(o). 

Lemna  2.2.1:  If  the  velocity  profile  satisfies  (2.2.27) 

and  for  some  n 2 0,  ^ erg  , 


x*(«i)  ^ x*(Qfa)  . (2.2.28) 

then  for  m 2 n, 

^ 2Vt  . (2.2.29) 


Proof:  We  shall  show  that 

X*(CTa)  - 2Vt  s X*(Oi)  s X*(Oa),  m 2 n.  (2.2.30) 
From  this  we  get,  by  (2.2.21)  and  (1.1.28), 

x*(Q(a)  - 2Vt  s x*(«i)  s x*(Oa)  + 2Vt  , 

which  Is  the  same  as  (2.2.29).  Equation  (2.2.30)  can  be  established 
for  m - n,  since 

X*(Q(a)  - X*(0fi)  - x*(Oa)  - *„(“!)  + T(v(Oi  )-v(Oa))  ^ 2Vt  , 

(2.2.31) 

by  (2.2.28)  and  (1.1.28).  We  suppose  that  (2.2.30)  Is  true  for  m, 
and  prove  It  for  m -f  I. 

One  half  of  (2.2.30)  Is  obvious: 
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X*  (Oi)  2 inf  x*(Qf) 

iotI  r Tit  n xn 

a€[U  .a^J 


(2.2.32b) 


i ®“P  4.  ®“P 

\a€[cifi,U  ] a€[u  .ofi  ] Qf€[Q'8,o'i] 


inf  x*(a) , sup  x*(of)  - inf  x*(a) 


a€[u;aa] 


oup  yji)  - A \Ji)\ 


s max  / 2Vt,  sup  x (a)  - inf 


a^[oa.c*il 


f x*(a)\ 

.c^i]  "*  > 


(2.2.33) 


x*(o)  - x*(p)  s 2Vt  for  a i p 


(2.2.34) 


Giv«n  <1  > 0,  Sa  > 0,  we  find  03  and  04  such  that 


03  € [oa.ofil,  x*(0f3)  ^ "up  x*(a)  - , 

oeCaa.cKi] 


«4  € [oa.ori],  x*(Qf4)  s inf  x (or)  + Ca  • 


(2.2.35) 


o^[aa.«i]  " 


For  the  other  half,  recall  from  the  algorithm  ( 1. 1. 12)-(1. 1. 14) 
that 


^ ®“P  + ’ 


(2.2.32a) 


X (or,  ) s sup  X (or)  - inf  x*(Qr) 

n^[aa>y^l  Q^[U  ,ori  J 
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If  QT3  a cx*.  by  (2.2.34)  and  (2.2.33)  we  get 

(“i)  ^ 2Vt  + ej  + 63  . (2.2.36) 

So  suppose  Qfa  < Q(^.  By  (2.2.30), 

X*(0f3)  - 2Vt  i;  X*(of3)  - 2vt  s x;J(ofi)  s X*(ct4)  ^ X*(03)  . 

(2.2.37) 

Using  (2.2.21)  and  (2.2.27),  we  find  that  this  leads  to 


**(0-3)  - 2Vt  s x*(o*) 


and  once  again  we  get  (2.2.36).  Since  e,  > 0 and  €3  > 0 are 
arbitrary,  we  get 


^ 2Vt 


(2.2.38) 


This  establishes  (2.2.30)  by  induction,  and  the  lemma  is  proven. 

Ihe  lemma  shows  that  if  X*(cxi ) and  X*(oa)  are  "close" 
then  ,c  will  x*(a.)  .„d  x‘(«.)  "for  „ i H 
this  as  giving  a Bound  on  the  width  of  numerical  shocks.  The  re- 
sult is  still  not  as  strong  as  we  desire,  since  we  have  no  idea 

“hocks  there  may  be.  One  may  envisage  the  algorithm 
(i. 1.12) -(1.1. 14)  as  unfolding  through  the  "interweaving"  of  vari- 
ous stre^.  It  then  appears  that,  once  two  such  streams  are  inter- 
woven, they  cannot  subsequently  become  unwoven.  Thus,  we  will 
naturally  want  to  study  families  of  such  interwoven  streams,  in  the 
hope  that  they  will  be  approximations  to  the  shocks  of  the  exact 
solution.  To  see  what  sorts  of  families  of  streams  we  should  con- 
sider, we  may  return  to  our  initial  heuristic  discussion  of  the 
exact  solution,  where,  we  recall,  we  searched  for  the  location  of 
sh^ks  from  "left"  to  "right",  mlnlmialna  expressions  like  (2.2.12) 
and  proceeding  downward  from  the  "top"  of  the  velocity  profile 
to  the  "bottom",  at  U“. 
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Accordingly,  we  might  consider  sets  of  the  form 


2 {oflof  € [U  a P ^ [U~,Qt]  such  that  x*(p)  - x*(u'*’)  s 0 


for  some  n s m - 1} 


(2.2.39) 


If  we  let 


9^  = inf  {v(Qf)|o(€B^} 


(2.2.40a) 


(2.2.40b) 


we  can  prove,  by  induction  in  m,  that 


- X^U-')  S T(v(u'^)-®  ) . 


(2.2.41) 


However,  it  may  be  that  a stream  becomes  interwoven  with  a 
second  stream  which  has  itself  become  interwoven  directly  with  the 
stream  tT*" . Such  streams  would  not  appear  directly  in  the  set  B„, 
introduced  in  (2.2.39).  A natural  extension  of  the  result  given 
in  the  last  paragraph  then  involves,  given  a set  DC  intro- 

ducing a set 


B^*(D)  = {o|Qf€  [U‘,U'*'1,  a p€  [u",or]  and  3 6€  D such  that  x^(B)-x^(6)  ^ O) , 


(2.2.42) 


We  observe  that 


=»  D 


(2.2.43a) 


B^(D)  - B^(D) 


(2.2.43b) 
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where  by  § we  mean  the  convex  hull  of  a set  S: 


S - {a|a  e € S and  a Y ^ S such  that  p s a s y}  . (2.2.44) 


Then,  we  can  prove  the  following:  If 


sup(x*(a)|a€  d)  - inf (x*(a) jaC  d}  £ d , (2.2.45a) 


then 


8upU*(a)lo'€  Bj^*(D)}  - inf[x*(a)l»€  Bj^*(D)] 


s max(d+T[sup[v(Q')lo'€  B *(D)-D]  - inf{v(o)la€  b^*(D)-d}], 


t[sup{v(cy)|q'€  Bj^*(D)}  - inf {v(o)1q’€  Bj^*(D)}}  } . (2.2.45b) 


We  are  getting  progressively  larger  sets  of  interwoven 
streams,  but  we  still  have  not  achieved  "closure"  of  these  fami- 
lies. The  next  natural  step  Is,  with  B„*(D)  defined  by  (2.2.42), 
to  generate  the  sequence  of  sets 


B<“> 


“ \*  (D)  , 

^o 


bC'-h-i) 


(B<™>) 


(2.2.46) 


Since  B^  ^ 3 B^  ^ and  B^  ^ C.  [u  ,U3^  m2  l,  we  get  conver- 
gence, and  can  define 


B * (D)  s U B^®^  . 
*0  m-i 


(2.2.47) 
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B *(B  *(D))  = By*(D) 


(2.2.48) 


With  this  much  as  background,  let  us  finally  define  the 
sets  which  will  appear  in  the  statement  of  our  principal  result. 
We  have  partitions  of  the  interval  [U",U+].  By  this,  we  mean 
that  there  are  ordered  sets  A„  ^ [U",U+]  of  points,  labeled  by 
the  parameter  or,  such  that  with  each  O € A„  there  is  associated 
a set  Dn(or)  with  the  following  properties: 


e s or  if  8 € D^(Qr)  , 


(2.2.49a) 


ori€A^,  as^A^,  > ofs  , 8i€D^(q),  8s  € D^(Os  ) =*  8,  > 8s  . 


We  further  require  that 


(2.2.49b) 


U'*’  € A V n 2 0 
n 


(2.2.50a) 


U D (a)  - [u‘,u''’]  V n 2 0 . 
o€An  " 


(2.2.50b) 


We  may  restate  (2.2.49b)  as 


Qfi  € A^,  oa  € A^,  ^ - <p 


(2.2.51) 


The  sets  A^  will  have  the  property  that 


Vi  ^ ^ Ao  - [u',u'*-l 


(2.2.52) 
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Thus , 

Dto(a)  - (a),  a e . (2.2.53) 


(2.2.53)  describes  the  partition  rig.  We  will  have  described  all 
the  partitions  if  we  can  show  how  to  get  from  to 

To  get  from  II  to  0 we  do  the  following:  We  form  the 

n n+i 

function  X*(o?)  in  accordance  with  the  algorithm  ( 1 . 1 . 12) - ( 1 . 1^14) 
aid  the  definition  (2.2.19).  Suppose  we  have  a set  DC  [U  ,U  ]. 
Given  D,  we  define  a number  ®jj(D)  ' 


a^(D)  = sup{q'|q'  G A^,  3 0 ^ D such  that  p € 0^(0-) } . 


(2.2.54) 

Now  we  define  another  set  by 

C^*(D)  s {ala€  [u',a^(D)l,  3 \€  A^,  3 0€  [u‘,y1,  and  3 6€  d such 

n 

(2.2.55) 

that  Of  € D (Y)  and  x*(0)  - x*(6)  ^ 0}  . 
n n n 


Note  that 


and 


C *(D)  =>  D 
*n 


where  we  use  the  notation  (2.2.44).  Also,  if 


» C^*(D)  , 

^n  n 


D (a)  n c„*  (D)  ^ <t>  , 

n A 

n 


D^(a)  c Cjj*  (D)  , 

n 


(2.2.56a) 


(2.2.56b) 


(2.2.57a) 


I 

I 

I 


then 
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and,  if 


^ 0s>  , 


then 


C^*  (Di)  C (Da)  . 

n n 

Next  we  generate  the  sequence 


(2.2.58a) 


(2.2.58b) 


r^'^D)  = C * (D)  , 

*n 

r^“^D)  - c * (r^^^D)).  •••  , 

^n 

r^“^'\D)  - (r^”>(D)),  •••  . 

n 

Since  r^“^(D)  C [u’.u'*’]  and  ^ T^"*^  V tn  a i. 

e * (D)  = u r^”\D)  . 

Xn  m-i 


(2.2.59) 


we  can  define 


(2.2.60) 


For  each  or  € A we  find 
n 


°n+i^“^  “ ^X*  • (2.2.61) 

n 


Clear ly^ 


S*  <°n+i<“>>  “ Vi  <">  • (2.2.62) 

n 

It  follows  from  (2.2^54),  (2.2.55),  (2.2.59),  and  (2.2.60)  that 

U ^ Dn+i  (l^)  and  U ^ *^n+i  ® ^ \ * 

subset  of  points  a € such  that  the  associated  sets  D (a)  are 

n n+i^  ' 

a disjoint  collection  of  sets  covering  [U",lr].  From  the  above,  we 

muat  have  IT*"  In  this  subset.  Tentatively,  we  call  this  subset  of 
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I 


points  A , . Since  we  shall  see  in  the  next  paragraph  that  A 

n+i  '^n+i 

thus  defined  is  determined  uniquely,  we  have  also  defined  the 
partition 

To  see  that  is  defined  uniquely,  we  proceed  by  induc- 

tion. It  is  clear  from  (2.2.52)  that  Aq  is  uniquely  determined. 
Assume  that  A^  is.  Now,  suppose  that  A^^^  is  not  defined  uniquely, 

so  that  there  are  two  sets,  A , and  k'  , satisfying  the  condi- 

’ n+i  n+i  o 

tions  of  the  definition.  Let  a ^ A , , a 4 A^ . . Since 

n+i  ^ n+i 

A . c A , a ^ A . Also,  it  is  clear  from  the  preceding  results 
n+i  n’n  ® 

that  a ^ iT*".  Then  € A such  that  a^^'  > a,  a 6 D . , ')i 

and  a^'^^  € a'  . Since  A'  A , a^^^  ^ A . If  € A 

n+i  . n+i  n’  n n+i 

and  a € (o'  O,  we  get  a ^ requirement  that  the 

sets  [d  (a)la€A  } be  disjoint.  Since  this  is  a contradiction, 

n+i  n+i  - . - ^ . . 

we  find  a^' ^ ^ A . . So  g a^^^  > with  a'^^  € A , C A and 

a'  € D (a(®)).  Since  a'  € D (a'^0  = Cy*  (D  (a'®^))  by 

n+ 1 n+i  * » n+T 

(2.2.62),  it  follows  from  (2.2.57)  that  D (a'^^)  c d . (a'^O. 

^ V n n+i 

From  (2.2.58)  we  get  Cy*  (D  (a'^')  c D (a^®^),  and  in  general 
/ \ ^ n+i 

r^'"^(D^(a<^))  c D^^^(a(=)).  Hence,  from  (2,2.60)  and  (2.2.61), 

D . (a(^))  D . (a^®^).  Since  a € D . (or^^^),  we  get  a € 

n+i  ' ' n+i'  ' n+i ' ® 

D .,(Qf^®7)  for  a € A , and  a'  ' € A , . This  gives  us  a contra- 
n+i  n+i  n+i  ® 

diction  with  the  required  disjointness  of  the  sets  {D^^^(a)la  ^ 

A ].  Accordingly,  A as  defined  must  be  unique. 

1*1  I n*ri 

The  sets  D^(a) , a € A^,  have  the  following  property:  If 

we  replace  i/’(x)  by 

i^(x)  E min(if  (x),Y),  V ^ a , 

then  for  0:Sms:  nX  (8)  for  8 ^ D (a)  is  unchanged. 

m n 
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Let  US  now  prove  the  following  lemma. 

Lemma  2.2.2;  Let  the  sets  An  and  the  sets  Dn(P),  P ^ An, 
be  constructed  according  to  the  prescription  above.  Then  for 
n 2 0, 


sup{x*(or)  lcif€  D^(3)}  - inf  (X^Cor)  (orC  D^(P)}  s TCsupfvCa)  lor€  Dj^O)] 


- inf  {vCcif)  lof€  D^(e)}] 


(2.2.64) 


Proof:  The  case  n » 0 is  trivial,  since,  by  (2.2.53), 


sup{x*(of)|or€  Eb(p)}  - inf  {x*(a)  la€  (p) ] = X*(p)  , (2.2.65a) 


«up(v(a)|of€  DbO)}  -inf{v(a)lo6  Db(P)J  = v(g)  . (2.2.65b) 


We  proceed  by  induction.  Thus  assume  (2.2.64)  true  for  n.  We  prove 
it  for  n + 1. 

As  in  the  proof  of  leoma  2.2.1,  we  have 


^ ®“P  • (2.2.66a) 

0€[o,U+l  " 


xL,  (“)  • (2.2.66b) 

P€[u-,a] 

By  the  construction  of  D (P)  with  p € A , If  a « p and  a ^ D . (P), 

is»  1 irT  ^ n+i 

then  x*(a)  > x*(V)  for  all  Y ^ (P)*  Likewise,  if  o 2 P and 

n n n*ri 

O 4 (P),  then  x*(of)  < x*(Y)  for  all  Y € D^^^^(P).  Hence  the 

induction  will  follow  for  (2.2.64)  if  we  can  show  that 
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suplx*(Qf)lo  G D^^(3))  - lnf{x*(a)lo  G D^^^O)} 


(2.2.67) 


s T[sup(v(a)la  G D (P)}  - lnf{v(a)la  G D (p)}]  . 

n+1 


First,  we  try  to  bound 


8up[x*(a)lo;G  (D^(P))}  - lnftx*(o)la  G c^*  (d  (p))] 
^ n 


(2.2.68) 


As  in  the  proof  of  lemna  2,2.1,  this  can  be  approximated  arbi- 
trarily well  by  x*(Pi)  - xS(Pa)  for  some  p^  G (D  (P))  and 

p2  € (D^(p)).  If  Pi  i Pa, 


x*(0i)  - x*(Pa) 


x;(Pi)  - x;(Pa)  + T(v(Pi)-v(Pa)) 


s (v(Pi)-v(Pa))  , 


(2.2.69) 


and  we  get  the  desired  bound. 

Suppose  Pi  < Pa.  We  can  write  Pi  G D^(Yi),  Pa  ^ D (Ya)t 
such  that  a Pi'  s Y^  and  pj^  s Ya  with  x*(P/)  s x*(6i),  x*(P^)  i 
Xjj(6a)  and  6i  G D^(P),  6g  G D^(P),  Clearly,  Pi  G (D^(P)) 
and  Pa  G C^*  (D^(p)).  Since  Pi  < Pa , Pi  C D^(P)  =»  p"  € D^(P). 

In  that  case. 


‘ ^ xJ(Pi)  - x;J(pa)  s 8uplx;(a)lcy  G D^(p)) 

(2.2.70) 

“ ^ 'r[8up{v(o)laGDj^(P)}  - inf {v(a)l aGDj^(p) } ] , 


on  use  of  (2.2.64). 


So  consider  the  case  Pa  G D (p).  Pi  ^ D„(p). 
* Pi » where  pi  Is  given  above.  Then  “ 


Suppose 
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- x*(3i')  s T(v(pi)-v(Pi)) 

(2.2.71) 

S T[8upU(a)la€Cj^*(D^(P))-D^(B))-lnf{v(a)la€c„*(D  (B))-D  (P)]]. 

Otherwise  3^  > Pj  . Since  P^  € Dj^CYi  ) , we  ^ have  Pj  € D^(Yi  ) . 

^ 8‘*p{x*(0f)|Q'€D^(YO3  - lnf{x*(a)|a€p^(Yi)),  Yi  < p. 

(2.2.72) 


‘*°"®  before,  we  approximate  the  right  hand  side  of 

(2.2.72)  by  ’ 


- 5{J(Oa)  | 

(2.2.73) 

+ T[v(Qf,  ) -v(Qf8 ) ] , QTi 6 D^(Yi  ) . jOae  D^(Yi  ) . 

When  Qr^  ^ erg , an  upper  bound  is 

Xjj(ai)-  X*(Qfa)  s T[sup{v(o)lor€D^(Yi)}  - inf  tv(a)  1o6d^^;Yi  ))]  , (2.2.74) 

i 

^ ®®”®  •*°‘*"‘*  directly  from 

(2.2.73).  In  either  case,  we  get  I 

*>i)  - x;j(Qra) 

s T[.up{v(o)|a€Cjj*(D^(P))-D^(B)}  - inf  {v(o)  |a€r*(D  (j3))-D  (P)}], 
n Aji  n ’ n 

since  D^(Yi)  c (D^(P))  and  D^(Yi ) H D^(p)  -0.  This  last 

bound  is  the  .ane"as  the  right  hand  aide  of  (2.2,71).  which  holds, 
therefore,  whether  s 3^  or  > 0i • T 


.•'i' 
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Since  Pa  ^ 8®t,  with  6^  given  as  above, 

x*(Pi)-  x;j(aa)-x*(Pi)-x;(Pi)+  x*(p,')-  x*(60+x*(60  -x*(Pa) 
s T[8up[v(cr)la€Cj^*(D^(P))-D^(P)}  - inf  (v(Qr)  lcif€Cj^*(D^(P))-D^(P))] 

+ t[sup{v(q')|o€d^(P)}  - inf  {v(Q')  la€D^(P)}] 
s T[sup(v(ar)la€Cj^*(D^(P))}  - inf{v(a)lo€c^*(D^(P))}]  . (2.2.75) 

Here  we  have  bounded  x*(6^ ) - x*(Pa)  in  the  same  way  we  bounded 
x*(Oi)  - x*(aa)  in  (2.2.73). 

We  have  treated  the  cases  P^  s pg  and  Pi  < Pa , Pa  ^ D^(P). 
Suppose  now  Pi  < Pa . Bs  $ case  ga  ^ Yi  » with 

Yi  defined  above.  Then  since  Pi  < Pg  and  Pi  € D^(Yi ) , we  have 
Bs  ^ above, 

x*(Bi)-  x*(Pa)«  T[8up[v(0')lcif€D^(Vi)}  - inf  {v(a)lo!€D^(Yi  )}1.  (2.2.76) 

Otherwise  Pg  > Yi  and  Pg  € D^(Ya).  Yi  < Yg  < P.  Choose 
Of  € D^(P)  such  that  v(or)  approximates  inf  [v(o')  jof^D  (P)}.  Then 
Of  > Pa  . So 

*J(«)-xJ(Ba)  - X*(of)  - X*(Pa)  + T[v(a)-v(Pa)]  s T[v(a)-v(Pa)]  . 

(2.2.77) 

With  3i  given  above,  we  have,  for  Pi  i Pi , 

x*(Bi)  - xJ(Bi')  « t[v(Pi)-v(Pi')1 

i T[eup{v(a)loi€D^(Yi)}  - inf(v(a)la€Cjj*(D^(P))}]  . (2.2.78) 
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For  ^ 3i  . bounded  just  as  we 

did  in  (2. 2, 72) -(2. 2. 74).  A suitable  bound  was  given  by  the 
right  hand  side  of  (2.2.74),  which  in  turn  can  be  bounded  by 
the  right  hand  side  of  (2.2.78).  Thus  (2.2.78)  holds  whether 
Pi  s Pi  or  Pi'  > Pi . 

Choosing  a 6 D^(P)  appropriately,  we  can  bound  **(of)  - 
x*(0s)  above  as  closely  as  we  desire  by  r[inf [v(a) |q€d^(P)}  - 
v(Ps)]«  Finally,  write 


x*(Pi)-  x*(Pa)”x^(Pi  )-x*(Pi  )+x*(Pi)-x*(6i  )+x*(6i  )-x*(Q')-Hc*(a)-x*(82) 
s T[sup{v(of)la€D  (Yi)l  - inf (v(a) lo€c„*(D^(P))} ] 

+ T[sup(v(cr)lci(€D^(P)}  - inf [v(of)  1o€d^(P)  ]] 

+ T[inf{v(a)la€D^(B)}  - v(Pa)] 

s T[sup{v(af)|a€Cjj*(D^(p))}  - inf  {v(a)  lo€Cj^*(D^(P) ) ]1  , 

n 

Here  we  have  used  the  fact  that  Pg  > Yi  , and  we  have  bounded 
Xn(6i)  - in  a now  familiar  way. 

Thus,  we  obtain 

sup{x*(a)lQf€Cjj*(Djj(p)))  - inflxJJ(Qf)la€c  (D^(p))} 
n ^n 

(2.2.79: 

s t[sup{v(o()1q(€Cjj*(D^(B))]  - inf{v(a)lo€c  *(D  (P))}]  . 

n *n  ” 

Now  we  turn  to 


sup{x*(o)la€r^“^(D^(P))}  - inf(x*(Qf)loi€r^“^(Djj(B))}  , 


(2.2.80) 
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where  F''  '(D^(P))  was  given  in  (2.2.59).  As  before,  this  can  be 
approximated  by 


x*(PO  - infU*(a)la€F^=>(D^(P))],  F^®>(D^(P))  . 

For  o e F^®\d^(P))  - r^^^(D^(p)),  y*(a)  > xj(6)  for  all  6 € D^(P). 
Thus , 

inf{x*(Qf)laer^^\D„(P))}-  inf[x*(cy)lo€r^'^D„(p))].  (2.2.81) 


Suppose  Pi  € r^^^(D^(P)).  Then  sup[x*(cr)  lo€r^^^ (D^(P) ) } can  be 
approximated  by  sup{x*(of)  |ot€r^'^  ^ (D^(p)) } . 

Consider  the  other  case,  p^  ^ F^^\d^(P)).  P^  ^ D (Yi). 
Then  g Pi  ^ Yi  such  that  x*(Pi)  s x*(6i  ) for  some  6j  € F^^^(D^(P)). 
Pi  ^ r<  = )(D^(p)).  If  Yi  e r(0(D^(e)),  «e  get  Pi  € r^'^D^(0)), 
which  is  a contradiction.  Thus  Yi  ^ F^^'(D^(P)).  On  the  other 
hand, 

pi€  r^'’^D^(P)),  Pi€  D^(Yi)  » Yi6  F^'’^D^(8)).  So 
Yi€  F^’*^(D^(P))  - r^^^D^(P)).  Similarly,  Pi  s Yi  , 

Yi^  r^®\o^(P))  - r^^^(D^(P)),  and  the  connectedness  of 
r^'^D^(P))  imply  that  p/  6 F<=^D^(P))  - F<'^D^(P)). 

Consider  the  case  3i  £ Pi  . 


xJ(Pl)  - xJ(Pi')  - X*(Pi)  - X*(Pi')  + T(v(Pi)-V(Pi'))  s T[v(Pi)-V(Pi')] 


ifi  TCsuptv(a)|o(€r^®>(D„(B))  - r^'^D^(P))} 

- lnflv(a)la€F<*>(D„(p))  - F^^^D„(P))}]  . 


(2.2.82) 
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In  the  other  case,  Since  Pj  € D (Yi ) we  have  Pi  € D (Yi ) . 

As  before,  we  get  in  this  case  " " 

»'*(Px)  " ^ t[sup(v(q')1o€d^(Yi)}  - lnf(v(Q')|Qf€D^(Yi )]]  , 

which  in  turn  is  bounded  by  the  right  hand  side  of  (2.2.82)  since 

Yi  € r^®>(D^(p))  - r<'\D^(P)). 

Thus, 

x*(Pi)  - inf(xJ(Qf)lcy€r<'^D^O))}  = x*(Pi ) - x*(Pi')  + xJ(Bi')  - x*(6i) 

+ x*(6i)  - inf[x*(a)lo(€r^^\D^(P))} 

s T[sup{v(a)|a€r^®\D^(P))-r^"\D^(B))}-inf{v(a)la€r(^^(p))-r^'^D^(P))}] 
+ T[sup{v(cy)lo(€r^'^(D^(P))}  - inf(v(a)lo€r^'\D^(P))}l 
s T[8up(v(Qr)la€r^^^D^(p)))  - inf  {v(cy)  la€r^®\D^(P)) } ] . (2.2.83) 


We  may  continue  in  this  way,  getting 


8up{x*(a)|a€r^“‘^(D  (p))3  - in£[x*(a)|o(€r^‘“^(D  (p))} 

(2.2.84) 

s T[gup{v(a)|a€r^”\D^(P))}  - inf  {v(o)  lcy€r^®\D^(P))  ]] 

and 

•up[x*(Qf)|a€Cj^*(Djj(P))3  - inftx*(cr)lo  €?^*(D^(P))J 

n n (2.2.85) 


s T[sup{v(a)|a€Cj^*(D^(P))}  - inf {v(a)|a  € Cjj*(D^(P))}] 


1 

1 

i 

1 

I 

I 

I 

] 

I 

1 

I 

I 

I 

1 

] 

I 


-78- 


THt  JOHNS  HOPKINS  university 

APPLIED  PHYSICS  LABORATORY 

MARYLANO 


This  establishes  (2.2.67)  by  induction 


and  the  lemma 


Thus,  for  V monotone,  t 
run  down  from  a peak  at  IT*"  to  a 
t|v(U+)-v(U')|  . A natural  conj 
The  total  width  of  all  shocks 
a minimum  U is  bounded  by 


lie  cocai  wiotn  ot  aiL  shocks  as  we 
. minimum  at  U”  is  no  more  than 
ecture  to  make  is  the  following: 
as  we  run  down  from  a maximum  IT*'  to 


where  the  variation  in  the  Interval  is  defined  by  (Reference  8) 


and  Tfis  a partition  of  as  follows 


We  leave  the  proof  (or  disproof)  of  this  conjecture  to  the  reader 
Since  a function  of  bounded  v.iation  is  differentiable  almost 
everj^ere,  it  is  clear  anyway  that  v of  bounded  variation  in 
[U  ,lr]  can  be  approximated  as  closely  as  we  desire  (in  the  maxi- 
mum norm)  by  a piecewise  constant  v.  And  for  a.u.c.  initial  data 
we  have  bounded  the  error  Induced  by  such  an  approximation  in 
lemma  2.1.3. 


£1.  situation  is  this,  then;  After  n time  steps  we  have  a 

profile  X (a)  and,  for  3 € A^,  shocks  spanning  the  range  of  veloci- 
ties In  the  Dn(0).  To  complete  our  analysis,  we  should  compare  our 
*2**2°*^"**^^  solutions  with  the  exact  solution  X (P)  given  in 
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Ui<Du“;n)  « Ui(Dif  ;C1)  + mr  J*  v(o)dtt,  0 s m s n , (2.2.87a> 

where 

0=  {(x,ar)lx2  X?(U'*’)  - nrv,  a€D^(0),  B € A^}  . (2.2.87b) 

Defining,  for  $ € A , 


(e)  = inf {ala€D^(8)}  (2.2.88) 

we  find  from  (1.1.24)  that 

-P  P 

nT  ^ v(Qr)do  » (X*^(Qf)-X*(o')  )do'.  (2.2.89) 

Using  the  result  in  leirana  2.2.2  that,  for  <y  € D (8) . B € A . 

n n 

sup{X*(a)|c«€D  (8)}  - T[sup{v(Qf)  |ck€d  (8)}  - inf  {v(o')  joCD  (8)}] 

‘•a*  n n 

i X*(or)  s inf(x*(Q')|o(€D^(8)J+T[aup{v(Q')lQ!€D^(8))-lnf{v(Qf)la€D^(8))]  , 


we  see  that  (2.2.88)  leads  to 


lx!(a)  - L x^(o)do  - -SrL  («)<tol 


e-D<  1)< 


0-D<  - 


(2.2.90) 


s T[sup{v(£y)|o(eD^(8))  - lnf{v(of)lQf€D^(8)}]  , 


for  all  8 € A„  such  that  D<  (8)  i-  8*  In  this  paragraph  and  the  fol- 
lowing ones,  one  may.  If  ha  desires,  consider  that  the  nonotonlcally 
non-increasing  function  (x)  has  bean  replaced  by  a histogram,  as 
described  in  (2.2. l)-(2.2.4) . In  that  case,  the  sets  A„  already 
introduced  and  Che  sets  Be  Co  be  incroduced  are  all  finiCe  and  1^(8) 

< 8 for  8 ^ An*  Our  final  error  escimace  will  be  independenc 
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pendent  of  how  closely  we  approximate  (x) 
will  hold  quite  generally.  Alternatively,  w 
ing  simple  bound  for  the  error  induced  in  th 
Xf.  and  X„  by  replacing  u°  (x)  by  the  histogra 
monotonically  non-increasing  and  satisfying 
( • »7).  For,  it  is  clear  in  this  case  that 


m6  S X < (m+l)6 


Given  the  exi 
a partition  9^  of  [u 
define  a set 


n (z.2.16),  we  can  effect 
For  each  Y ^ we 


a set  ^ [U  .U  ] such  that  for  \ ^ Bj. 
disjoint  covering  of  [u  '.U'*’].  We  let 
Y}.  Just  as  we  proved  that  the  sets  A_ 
defined  uniquely,  we  can  show  that  the 
rained.  As  we  defined  D<-(9)  in  (2.2.88) 
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The  analogue  of  (2.2.87)  in  this  case  is 


Ui  (Du(.,t');n)  =»  Ui  (Du(-,0);Q)  +t'  r v(a)da,  Ost'st, 

Et(V) 


where 


(2.2.95a) 

(12  {(x,or)l  X 2 Xo(u"^)  - tv,  at  E^(y) , y€  . (2.2.95b) 


Lemma  2.2.3;  If  u° (x)  - u(x,0), 


A C B 
n nT 


(2.2.96) 


Proof:  It  follows  from  (2.2.10)  and  (2.2.19)  that  X^(a)  = 

Xo(a).  The  lemma  is  trivially  true  for  n = 0.  We  will  prove  the 
lemma  for  n i 1 by  contradiction.  So  suppose  P € A„  and  8 € E (Y) , 
Y > P,  Y ^ First,  consider  the  case  E^(Y)  < 8-  Over  the’^'Set 


= {(x.of)l  X 2 Xo(U'*')  - nxv,  E^(y)  < a < 8}  , (2.2.97) 


we  can  have  no  gain  of  Ui  (fl)  in  the  approximate  solution,  and  we 
can  have  no  loss  of  (fl)  in  the  exact  solution,  beyond  that  which 
flows  across  the  boundary  at  x - Xo(ir*")  - nrV.  Thus, 


B 

(B-E^(Y)X  (Y)  2 s [Xo(cy)+ nTv(a)]do  , (2.2.98a) 


E<(Y) 


3 e 

**  X*(a)do(  s j"  [Xo(o:)  + nrv(a)lda  , 


E<(Y) 


E<(Y) 


and 


x;(0  ) i X^^(Y)  . 


Similarly,  over  the  set 


(2.2.98b) 


(2.2.98c) 


1 
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s {(x,a)l  X ^ Xo(u'*')  - nxv,  p < a < y]  , (2.2.99) 

we  can  have  no  loss  of  (Q)  in  the  approximate  solution  nor  gain 
of  (Q)  in  the  exact  solution,  except  for  the  flows  across  x = 
Xo(U+)  - nTV.  This  gives 


X*(a)do'  2 [Xo(a)+  nTv(a)]dQ;  , 

0 0 

(2.2.100a) 

Y 

(Y-8)X  (Y)  ^ J"  [Xo(o')  + nTv(a)]do'  , 

nr  g 

(2.2.100b) 

X*(P^)  ^ X^^(Y)  . 

(2.2.100c) 

Since  monotonicity  cf  X*  requires  X*(P  ) s X*(P  ),  the  inequalities 
in  (2.2.98)  and  (2.2.100)  must  be  equalities  and  X*(P)  = XnT-(Y). 

In  particular,  from  (2,2. 100a, b)  and  the  monotonicity  of  X^  we  con- 
clude that 


I 

I 

I 

I 

I 

I 


But  this  violates  the  definition  of  0 € A , and  we  get  a contra- 
diction. " 

Next,  consider  the  case  E<(Y)  = 0.  The  results  (2.2.100) 
still  hold.  On  the  other  hand,  from  (2.2.16)  and  our  construc- 
tion ( 1. 1. 12) -( 1. 1. 14) , we  still  must  have 


XnxCY)  - X^^(0)  2 Xo(8)  + nrv(0)  = X^(0)  + nTv(0)  2 X*(P) 

(2.2.101) 

if  0 € Ajj,  the  inequalities  in  (2.2.100)  must  become  equalities, 
and  we  get  a contradiction  as  before. 

Thus,  if  0 € Ajj,  we  also  have  0 € . Two  cases  may 

arise:  (i)  Dn(B)  - E„t(P);  or  (11)  3 y € 0^(0)  with  V < 0 and 
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Y ^ ®nT’  first  case,  the  same  considerations  as  above  show 

tha  t 

P 

J X*(o;)dcY  = (P-D  (P))  X (P)  if  D (p)  < p , (2.2.102a) 

D<(P)  " < nr  < 

X*(P)  = X„^(P)  if  D<(P)  =P  . (2.2.102b) 


Either  way, 

inf{x*(a)lQf€  D^(P)}  ^ X^^(P)  s sup{x*(a)  ja  € D^(P) } . 

(2.2.103) 

In  the  second  case,  if  D<;(P)  < Y P>  our  arguments  based 
on  the  flow  of  indicate  that 


Yr  i*P 

J X^(Q')do'  £ ; X^^(0')da  (2.2.104a) 


and 

r Y 

J X*(Qf)dQ'  s I X (a)  da  . (2.2.104b) 

D<(P)  D<(P) 


From  (2.2.104a)  and  the  monotonicity  of  both  X*  and  X , we  conclude 
/ n nT 


lnf{x*(a)laE  D^(p)}  s X^^(Y)  . (2.2.105a) 


Likewise,  (2.2,104b)  leads  to 


X„t<Y)  s suplx*(a)la6  D^(P)]  . (2.2.105b) 


Finally,  suppose  Y - D<(P)  < p.  The  result  (2.2.105a)  still  holds. 
On  the  other  hand,  a direct  check  of  (2.2.16)  and  the  construction 
(1. 1. 12)-(1. 1. 14)  shows  that 
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X„^(Y)  - Xo(Y)  + nrv(Y)  “ X?(Y)  + mtv(Y)  5 X*(Y)  . 

(2.2.106) 

In  all  cases  considered,  we  have  shown  that,  for  Y ^ D^(P) , 

lnf{x*(of)|cY€  Dn  (p)}  s X^(Y)  s sup{x*(a)|a€  D^(P)}  . 

(2.2.107) 

Combining  (2.2.107)  with  (2.2.64),  see  that  we  have  proven  the 
following  theorem. 

Theorem:  For  Y ^ ®n(P)» 

1 X*(Y) -X^^ (Y) I ^ T[sup{v(Qf)la€  d^(P)3  - inf (v(a) la€  D^(P)] ] . 

(2.2.108) 
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